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BRIEF  OUTLINE  OF  RESEARCH  FINDINGS 


This  final  report  summarizes  the  work  on  ARO  Grant  Number 
DAAG29-79-C-0082  from  1  May  1979  to  30  October  1982.  The  purpose  of 
this  grant  is  to  develop  theory  and  techniques  for  small  antennas 
mounted  on  structures,  for  microstrip  antennas,  and  for  k-pulse 
applications. 

Much  of  the  research  effort  during  this  period  is  described  in  the 
published  papers  which  are  reproduced  in  the  Appendices. 

First  let  us  summarize  our  accomplishments  related  to  small 
antennas  mounted  on  structures. 

The  theory  and  computer  programs  were  developed  for  plane-wave 
scattering  by  an  infinitely  long  elliptic  cylinder  with  arbitrary 
surface  impedance.  The  solution  has  the  form  of  a  series  of  Mathieu 
functions.  The  calculated  data  show  excellent  agreement  with  GTD 
calculations. 

We  investigated  the  characteristic  modes  of  Garbacz  and  Harrington 
for  conducting  strips,  strip  gratings,  thick  tubular  dipoles,  and 
periodic  collinear  arrays  of  thick  tubular  dipoles.  For  the  conducting 
strip,  it  was  found  the  characteristic  modes  coincide  with  the 
eigenmodes  of  the  degenerate  elliptic  cylinder  expressed  in  terms  of 
Mathieu  functions.  For  periodic  structures  such  as  the  strip  grating 
and  the  collinear  array  of  dipoles,  when  the  spacing  is  less  than  the 
wavelength,  a  single  characteristic  mode  is  induced  by  a  plane  wave  with 
normal  incidence. 


In  most  applications  It  Is  difficult  to  obtain  more  than  three  or 
four  of  the  characteristic  modes.  In  the  complex  frequency  domain,  new 
techniques  are  required  for  characteristic  modes  since  the  Impedance 
matrix  Is  no  longer  positive  definite.  For  these  reasons,  we 
investigated  the  conjugate  modes  of  Inagakl.  Even  in  the  complex 
frequency  domain.  It  Is  always  easy  to  obtain  a  complete  set  of 
conjugate  modes.  We  found  that  the  conjugate  modes  reduce  to  the 
natural  modes  at  the  poles  of  the  conducting  body.  For  the  strip  and 
the  tubular  dipole,  graphs  were  plotted  to  Illustrate  the  lowest-order 
characteristic  modes  and  the  conjugate  modes. 

We  applied  the  double  Fourier  transform  to  develop  the  theory  of 
thin-wire  and  thin-strip  antennas  located  on  or  near  the  air-earth 
Interface.  Previous  Investigators  have  used  the  Fourler-Bessel 
transform  for  these  problems,  but  the  double  Fourier  transform  has 
significant  advantages.  For  a  thin-strip  antenna  near  the  air-earth 
Interface,  our  technique  requires  double  numerical  Integration  whereas 
the  earlier  methods  require  five-fold  numerical  integration.  The 
Fourler-Bessel  method  encounters  a  singularity  when  the  antenna  is  on 
the  Interface,  but  the  double  Fourier  transform  technique  1s  well 
behaved  In  this  situation. 

Next,  we  applied  the  double  Fourier  transform  to  calculate  the 
current  distribution  and  Impedance  of  a  horizontal  strip  dipole  as  a 
function  of  distance  from  the  air-earth  Interface.  The  calculations 
Include  the  cases  where  the  dipole  Is  on  the  Interface,  burled  In  the 
earth,  and  elevated  above  the  Interface.  Although  the  Impedance  varies 


rapidly  with  distance  when  the  dipole  is  near  the  Interface,  it  appears 
that  the  impedance  is  continuous  across  the  Interface. 

We  also  considered  a  skewed  strip  dipole  near  the  air-earth 
interface.  The  skewed  dipole  is  neither  vertical  nor  horizontal,  but  Is 
oriented  at  an  arbitrary  angle  with  respect  to  the  Interface.  The 
double  Fourier  transform  proved  quite  advantageous  for  this  application. 
With  this  approach,  the  number  of  Sommerfeld-integral  evaluations  Is 
reduced  by  a  factor  of  ten  in  comparison  with  the  conventional  Fourler- 
Bessel  transformation. 

Next  let  us  summarize  our  research  accomplishments  related  to 
microstrip  antennas. 

We  applied  the  double  Fourier  transform  to  develop  the  theory  of 
microstrip  antennas.  In  calculating  the  self  Impedance  of  a  rectangular 
microstrip  antenna,  our  new  computer  code  reduces  the  computational 
expense  by  a  factor  of  ten. 

In  many  applications,  the  patch  model  of  a  microstrip  antenna  will 
Involve  trapezoidal  patches.  Therefore,  we  considered  the  mutual 
Impedance  of  two  trapezoidal  patches  and  developed  the  theory  far  enough 
to  show  that  the  transform  approach  Is  quite  feasible  for  this 
application. 

Finally,  let  us  summarize  our  research  related  to  k-pulse 
applications.  We  extended  the  k-pulse  theory  to  Include  the  effects  of 
the  generator  Impedance  at  the  transmitting  antenna  and  the  load 
Impedance  at  the  receiving  antenna.  For  numerical  calculations,  we 
considered  two  center-fed  cylindrical -wire  antennas  with  the  following 


generator  and  load  Impedances:  0,  50,  916,  and  •  ohms.  We  used  the 
frequency -domain  moment  method  to  calculate  the  antenna  Impedance  and 
effective  length  at  many  frequencies.  The  effective  length  was 
calculated  for  broadside  incidence  and  also  for  an  angle  of  45  degrees 
away  from  broadside.  From  these  data  we  employed  the  Fourier  transform 
to  calculate  the  following  k-pulse  waveforms:  the  generator  voltage 
Vg{t),  the  incident  far-zone  field  waveform  Ei(t),  and  the  received 
voltage  waveform  Vp(t).  A  manuscript  is  being  prepared  to  present  these 
new  results. 
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Mutual  impedance  between  vertical  dipoles  over  a  flat  earth 
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The  fields  of  a  vertical  dipole  and  the  mutual  impedance  between  vertical  dipoles  over  a  flat 
earth  are  expressed  in  terms  of  a  Sommerfeld  integral. 


1.  INTRODUCTION 

The  use  of  vertical  wire  antennas,  or  arrays  of 
vertical  wires,  over  the  earth  is  common.  This 
problem  can  be  analyzed  by  using  the  method  of 
moments  [Harrington,  1968] ,  where  the  effects  of 
the  flat  earth  are  accounted  for  by  modifying  each 
element  in  the  impedance  matrix.  The  computation 
of  the  modification  to  a  single  element  in  the  matrix 
involves  the  numerical  evaluation  of  a  Sommerfeld 
integral  [Sommerfeld,  1964] .  With  sinusoidal  bases 
and  point  matching,  it  has  been  shown  [Miller  el 
al.,  1972)  that  a  single  evaluation  of  the  Sommerfeld 
integral  suffices  for  each  impedance  element.  Miller 
et  al.  report  a  computation  time  of  20  s  on  a 
CDC-6600  computer  for  impedance  of  a  vertical 
half-wave  dipole  over  the  flat  earth.  In  view  of 
this  computational  expense,  they  present  data  based 
on  a  S-segment  antenna  model,  although  these  data 
differ  by  8%  from  those  with  a  21-segment  model. 
They  present  the  theory  based  on  enforcement  of 
the  integral  equation  on  the  axis  of  the  wire. 

In  this  paper  we  present  the  general  formulation 
and  show  that  a  single  evaluation  of  the  Sommerfeld 
integral  suffices  for  almost  any  choice  of  basis  and 
testing  functions.  Instead  of  matching  on  the  axis, 
we  present  the  theory  appropriate  for  enforcement 
of  the  integral  equation  on  the  surface  of  the  wire. 
More  accurate  data  are  presented  for  the  impedance 
of  a  center-fed  vertical  wire  versus  height  above 
the  earth,  and  it  is  found  that  our  calculations  are 
significantly  faster  than  those  quoted  above. 

2.  FIELD  OF  VERTICAL  DIPOLE 

For  an  electric  line  source  located  on  the  z  axis 
in  free  space  the  vector  potential  has  only  a  z 
component  given  by 

Copyniiit  Cl  ltT9  by  lb*  AiscncM  0«opbytical  Uruo*. 


A  = 


f  '5 


dz' 


(1) 


where  z  \  and  z\  denote  the  end  points  of  the  line 
source  and  /(z')  represents  the  current  distribution 
with  the  time  dependence  e^“'  suppressed. 

By  expanding  the  spherical  wave  exp  (-jkR)/R 
in  a  spectrum  of  cylindrical  waves,  we  obtain 


f"  F(y)e-” 

pAtpp) - dp 

(2) 

Jo 

y'  =  p'  -  k^ 

(3) 

F(7)=  j  '  l(z')e-'‘  dz' 

(4) 

where  z  is  less  than  zj.  The  electric  field  intensity 
is  given  by 


£.  =  CJ  gV,(ppXf/y)e-"dp  (6) 

where  C  =  —jy\/(47tk)  and  n  is  the  impedance 
of  free  space. 

Now  let  the  xy  plane  be  the  air-earth  interface. 
With  the  field  in  (6)  incident  on  the  interface  the 
reflected  field  is  given  by 

£'.  =  C  j  p’  /ty„(P,,K£/-Y)e  •"  dp  (7) 

where  the  reflection  coefficient  £  is  a  function  of 
1  and  (he  parameters  of  the  earth;  (p,  4i,  z)  denote 
cylindrical  coordinates  with  origin  on  the  interface 
and  line  source  on  the  z  axis. 

Now  consider  a  tubular  vertical  dipole  with  radius 
‘a,’  current  distribution  t(z'),  and  surface  current 
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density  uniformly  distributed  around  the  tube.  The 
reflected  field  fur  this  source  is  found  from  (7) 
by  using  the  addition  theorem  for  the  Bessel  function 
and  integrating  around  the  surface  of  the  tube: 

£:  =  cj  (8) 

where  p  is  measured  from  the  axis  of  the  tubular 
source. 

3.  MUTUAL  IMPEDANCE  OF  VERTICAL  DIPOLES 

The  mutual  impedance  between  tubular  vertical 
dipoles  is  given  by 

/,  and  /j  denote  the  terminal  currents,  and  the 
integration  extends  over  the  cylindrical  surface  of 
dipole  2.  Dipole  2  has  radius  a, ,  current  distribution 
/i(z),  and  uniform  surface  current  density  7}  around 
the  tube.  The  z  component  of  the  field  of  dipole 
1  is  given  by 

£.  =  £?+£;  (10) 

where  the  first  term  denotes  the  free-space  field 
and  the  second  term  is  the  field  reflected  from 
the  interface.  From  (9)  and  (10), 

=  +  (II) 

where  the  first  term  is  the  mutual  impedance  in 
free  space  and  the  second  term  is  the  change  in 
mutual  impedance  arising  from  interface  reflection. 

Using  (8)  for  the  reflected  field,  invoking  the 
addition  theorem,  and  integrating  around  the  surface 
of  dipole  2,  we  find 

=  (  3’(«/7) 

Jo 

•  7o(Pfl,)y„(Pa,)y„(P</)£C</p  (12) 

C(7)=|  l,(z)e''‘dz  (13) 

The  horizontal  distance  between  the  dipole  axes 
is  denoted  by  d,  and  z,  and  z,  are  the  end  points 
of  dipole  2. 

In  moment  method  applications  the  current  dis¬ 
tributions  and  /„(z)  represent  the  ba.si.s  func¬ 
tions.  and  (12)  determines  Z'„^.  For  any  simple 
basis  functions  it  is  apparent  from  (4)  and  (13)  that 


the  functions  F  and  G  can  be  obtained  in  simple 
closed  form.  Thus  we  see  from  (12)  that  a  single 
evaluation  of  the  Sommerfeld  integral  is  sufficient 
to  calculate  each  of  the  mutual  impedance  terms 
In  the  sinusoidal  Galerkin  formulation  the 
basis  and  testing  functions  are 


/-(r') 


sink  (A,  -  |z'  -  z;;.|) 
sin  kh^ 


(14) 


l.(i)  = 


sinlt(*.  -  |z  -  z.l) 
sin  kh. 


(15) 


where  and  denote  the  segment  lengths.  Func¬ 
tion  m  extends  over  the  two  segments  which  inter¬ 
sect  at  z'^,  and  function  n  extends  over  the  two 
segments  which  intersect  at  z. . 


4.  NUMERICAL  RESULTS 

Figure  1  illustrates  the  impedance  of  a  half-wave 
center-fed  vertical  wire.  Here  h  denotes  the  distance 
from  the  air-earth  interface  to  the  center  of  the 
wire.  These  data  were  calculated  with  the  sinusoidal 
Galerkin  formulation.  The  wire  was  divided  into 
30  segments,  so  there  were  29  equations  and  29 
unknowns.  In  Figure  1  the  trends  are  similar  to 


Fi(.  I.  Impedance  of  cenier-retl  vertical  wire  versus  height 
above  ibe  earth. 
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those  in  Figure  2  of  Miller  el  al.  [1972],  but  the 
dau  agree  more  closely  with  their  2 1  -segment  results 
rather  than  the  S-segment  results  they  present. 

If  the  20-s  computation  time  quoted  by  Miller 
et  al.  applies  to  their  S-segment  model,  our  corres¬ 
ponding  6-segment  calculations  are  100  times  faster. 
If  it  applies  to  their  21-segment  model,  our  22-seg¬ 
ment  calculations  are  12  times  faster.  We  employed 
a  simple  Newton-Cotes  integration  as  opposed  to 
their  adaptive  integration  for  the  Sommerfeld  inte¬ 
gral.  In  addition,  our  free-space  impedance  terms 
are  in  simple  closed  form,  whereas  theirs  may 
require  numerical  integration. 

Surface  matching  is  slower  than  axis  matching, 
and  the  two  techniques  give  essentially  the  same 
results  for  wires  with  small  radius  as  in  Figure  I. 
For  larger  radii  such  as  the  often  quoted  results 


for  a  =  0.007  X,  however,  surface  matching  is 
required  for  accurate  results  and,  indeed,  for  con¬ 
vergence. 
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On  the  Edge  Mode  in  the  Theory  of  TM  Scattering  by  a  Strip  or  Strip  Grating 

JACK  H.  KICHMOND.  iriiow.  iin, 


Abstraci — Considrr  a  plan*  wave  incideni  on  a  perfectly  conduciinx  our  calculations  with  three  unknowns  show  excellent  agree- 


slrip  (or  strip  ,  i  •a),  and  let  the  incident  electric  vecliir  be  parallel 
with  the  edKes  ui  the  strip.  If  the  edite  mode  is  included  among  the 
basis  functions,  it  is  found  this  greatly  improves  the  convergence  of 
lh«  moment-method  solution.  Numerical  data  are  included  for  the 
rellecllon  coefficient  of  (he  strip  grating.  To  correct  an  error  in  the 
previous  literature,  the  rigorous  solution  is  tabulated  for  broadside 
bacfcscallcr  from  a  single  strip. 

1.  INTRODUCTION 

IBLIOGRAPIIIES  on  strip  gratings  arc  available  in  1 1 1 
and  [21,  and  a  chapter  on  scattering  by  a  strip  is  included 
in  (3).  For  a  related  geometry.  Minor  and  Bollc  [4|  include 
the  edge  modes  in  an  efficient  formulation  for  shielded  micro¬ 
strip.  in  scattering  by  a  strip,  Shafai  and  El-Moazzen  |S|  use 
a  transformation  to  treat  the  edge  singularities  in  the  current 
distribution.  Tsai  et  al.  [6|  sample  the  current  function  at 
progressively  smaller  intervals  near  the  edge.  Wilton  and 
Govind  [7|  incorporate  the  edge  conditions  with  subsectional 
bases.  In  (61  and  (7|,  the  objective  is  to  obtain  more  accurate 
calculations  of  the  current  distribution.  Their  techniques 
evidently  did  not  improve  the  convergence  of  the  solution. 

In  this  paper  we  develop  the  theory  of  scattering  for  a 
strip  and  a  strip  grating.  When  the  basis  functions  include  the 
edge  mode,  it  is  shown  this  greatly  improves  the  convergence 
of  the  moment-method  solution  with  entire-domain  expansion 
functions.  For  a  strip  with  width  equal  to  the  wavelength. 

Manu.scripl  received  Nnvenibet  26,  1979;  revised  May  27,  1980 
This  work  was  supported  in  part  by  Contract  l>AAG29-79<'4J0K2  be¬ 
tween  Oepartrnent  of  the  Army.  II.S.  Army  Kescarcli  Ollicc,  and  the 
Ohio  Slate  University  Research  i  nundation. 

The  author  is  with  the  KiectroScicnce  l.abutatnry,  the  Ohio  State 
University, Columbus,  OH  4)212. 


inent  with  the  current  distribution  obtained  in  |7|  with  2$ 
unknowns.  For  a  strip  grating  with  normal  incidence  and 
spacing  less  than  the  wavelength,  it  is  found  the  current 
distribution  has  the  constant  phase  property  of  a  character¬ 
istic  mode  ( 8 1 ,  ( 9 1 . 

Although  our  technique  shows  some  similarity  to  the 
“asymptotic  anticipation  method"  of  Neureuther  and  Zaki 
(I.S|,  we  introduce  several  improvements.  In  their  discus¬ 
sion  of  plane  wave  scattering  by  a  wire  mesh.  Hill  and  Wait 
( 16)  use  related  techniques  to  improve  the  convergence. 

We  consider  a  grating  or  a  strip  with  perfect  conductivity, 
infinite  length,  and  innnitesimal  thickness.  The  time  de¬ 
pendence  is  assumed  and  suppressed  At  first  the  cur¬ 

rent  density  on  the  strip  is  expanded  in  a  Fourier  series,  and 
a  system  of  simultaneous  linear  equations  is  developed  for  the 
Fourier  coefficients.  Then  the  current  distribution  is  expressed 
as  the  sum  of  an  edge  mode  plus  a  Fourier  series,  and  it  is 
found  that  the  simultaneous  linear  equations  can  be  obtained 
readily  from  the  previous  expressions.  The  next  section 
considers  a  strip  grating  with  oblique  incidence,  but  for  the 
sake  of  brevity  (he  remainder  of  the  paper  considers  normal 
incidence 

II  TUI  ORY  OF  STRIP  GRATING  WITH  OBLIQUE 
INGIDI  Nfi; 

Consider  a  time-harmonic  plane  wave  in  free  space  to  have 
oMique  incidence  on  a  periodic  planar  array  of  thin  con¬ 
ducting  strips  as  in  Fig.  I .  liach  strip  has  width  w  -  2h  and  r 
denotes  the  spacing.  The  ;  axis  coincides  with  the  axis  of 
one  of  the  strips,  and  the  incident  electric  field  intensity  is 

=  (1) 
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of  simultaneous  linear  equations: 

/  =  0.±1.±2,  -  .  (12) 

/ 

ikv 


^il  ~  -  ^  ^in  *^//i/T/i 

2rM'  . 


(13) 


Fig.  1.  Plane  wave  with  oblique  incidence  on  strip  grating, 
where  k  =  uy/Jie.  The  scattered  field  is 


where  V/  =  0  except  Kq  -  Eq- 

The  reflection  and  transmission  coefficients  of  the  strip 
grating  are  given  by 


(2) 


R  =  Of) 
T=\+R. 


(14) 

(15) 


jun 

2  =  /'fr  A  X  _ yi 


(3) 


To  obtain  numerical  data,  we  truncate  the  series  in  (12)  and 
use  matrix  inversion  to  solve  these  equations  for  the  constants 
b/.  In  view  of  the  basis  and  testing  functions  used,  (12)  and 
(4)  (13)  may  be  considered  the  “Fourier-Galerkin  formulation." 

We  usually  expect  the  Galerkin  formulation  to  converge 
faster  than  a  point-matching  solution.  In  this  application. 
Let  us  expand  the  surface-current  density  on  the  strips  however,  the  two  techniques  show  equal  performance.  The 

only  advantage  of  the  Galerkin  method  is  that  it  yields  a  sym¬ 
metric  matrix  Z,/.  It  is,  however,  a  distinct  advantage  to  treat 
the  constants  b;  as  the  unknowns  in  (12),  rather  than  using 
the  On  as  the  unknowns  as  in  (IS). 


Tn‘  =  (*  tin  0  +  2/nr/s)*  —  k 

where  sgn  (x)  =  -  I  for  jr  <  0,  and  -H  for  x  >  0. 

Let  us  expand 
in  a  Fourier  series: 

From  the  boundary  condition  on  tangential  //, 

-/O.SkfjJb/e'^'’,  -h<y<h 
0,  in  the  apertures 


(5) 


111.  STRIP  GRATING  WITH  NORMAL  INCIDENCE 

For  normal  incidence  (0  =  0),  it  is  convenient  to  express 
the  scattered  field  as 


(6) 


ff*  =  zEa  2  cos  (2n»r>>/s) 


(16) 


where  rj  =  y/ti/f.  We  multiply  both  sides  of  (6)  by  e  and  the  current  density  as 

and  integrate  over  the  region  -s/2  <>  <s/2  to  obtain 


tffl  =  -»7  wbo/(2s£'g  cos  0) 
-jkV 


(7) 


y  =  i  2  cos  (litylh). 
0 


(17) 


On  = 


2rfnRo 

h 


'2 


(8)  By  enforcing  the  boundary  condition  on  tangential  //  (with 
testing  function  cos  2n7ty/s).  we  find 


Fi, 


.^=  I  g/ix(l/w-n/s)y  jy 

J-h 


(9) 


-/kt1>ve„  ^ 

<*n  -  ~  Zj  °filn 

Onto  0 


(18) 


By  forcing  tangential  E  to  vanish  on  the  conducting  strips,  we 
obtain  from  (1)  and  (2) 

=  -h<y<H.  (,0) 

We  multiply  both  sides  of  (10)  by  t-R’’*!''  and  integrate 
over  the  region  -h  <y  <h  to  obtain 


f/i/i^CI/'v)  /  cos  (Iiry//i)  cos  (2ntry/s)  dy  (19) 

A) 

where  Vn  =  2  except  Cq  =  1 .  By  forcing  tangential  f  to  vanish 
on  the  strips  we  find 


2®n  cos(2/nry/s)  =  — I,  ~h<y<h. 


(20) 


-w, 

0, 


lf/  =  0 
if  1*0 


(II)  Multiplying  (20)  by  cos  iny/h  and  integrating  over  the  range 
0  <y  <b  yields 


where  an  asterisk  denotes  the  complex  conjugate  and  Ff„  is  •• 
defined  by  (9).  In  (II)  let  us  replace  a„  with  the  equivalent  2/  = 

quantity  from  (7)  and  (8).  This  yields  the  following  system  n<u 


-0.5, 

0, 


if  /  =  0 
ifi*0. 


(21) 
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From  (18)  and  (21), 


N 


1  =  0, 1,2.  •./V 

1-0 

(22) 

m 

Zii  =  (2ikriw/s)  2  enGi„G,„iy„ 

(23) 

n=0 


where  Vj  =  0  except  Fo  =  A'q.  The  impedance  matrix  2,/  is 
symmetric.  The  reflection  coefficient  is  given  by  fi  =  Cq,  and 
do  and  bo  are  related  by  (7). 

If  the  incident  field  is 

i'  =  cos  (2miry/f), 

the  solution  is  given  by  (16)-(19)  and  (22)-(23)  with 

Vj  =  2EoGi„. 

IV.  STRIP  GRATING  WITH  EDGE  MODE 

The  current  distribution  J(y)  has  singularities  at  the  edges 
of  a  strip.  Thus,  when  J  is  expressed  as  a  Fourier  series,  the 
series  has  slow  convergence  and  one  must  invert  a  rather 
large  matrix  to  obtain  an  accurate  solution.  To  improve  the 
situation,  we  may  express  J  as  the  sum  of  an  edge  mode  (to 
take  care  of  the  singularities)  and  a  Fourier  series.  In  this 
case  the  Fourier  series  should  converge  rapidly  since  it  repre¬ 
sents  a  well-behaved  function. 

For  normal  incidence,  a  suitable  edge  mode  for  the  current 
distribution  on  a  strip  is 


2 


(24) 

Je  =  (2e/w)  2  COS  Q^y/h) 

(25) 

0 

f" 

^e=  /  ■fe(y)tly 

k 

(26) 

'-/i 


where  Joiln)  denotes  the  Bessel  function.  The  function  7,  has 
the  correct  behavior  in  the  neighborhood  of  each  edge  of  the 
strip. 

For  normal  incidence,  let  us  express  the  current  density 
on  the  strip  as  the  sum  of  the  edge  mode  and  a  Fourier  series: 

N-  I 

2  O  cos  (/try /A).  (27) 

0 

Using (2S)  for/,  and  comparing/ 1  7)  and  (27),  we  find 

bl  =  c,  +  (Ie/w)e,Jo(ln).  (28) 

From  (22)  and  (28), 

c/v  2^  2  c,7.i,=  (29) 

/•o  /-o 

where  <  ff  =  /,/w  denotes  the  amplitude  of  the  edge  mode. 
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Equation  (29)  can  be  rearranged  as  follows: 

N 

'^c,Zii'=Vi,  1  =  0,  1,  2,  • 

1=0 

(30) 

Zg/  =  z,/ ,  for  /  <  yv 

(31) 

nt-0 

(32) 

where  2<|  and  2,„,  are  given  by  (23)  and  Vi 

=  0  except  Kq  = 

Eq.  The  matrix  2,7'  is  not  symmetric. 

We  could  obtain  a  symmetric  matrix  by  using  the  edge 
mode  (instead  of  the  last  cosine  term)  as  a  testing  function 
with  (20)  to  derive  the  last  equation  (i  =  N)  in  (30).  In  this 
ease  2;v„'  =  and  K/  =  0  except  Vq  ~  ~  *0 

ever,  slightly  better  convergence  is  obtained  with  the  unsym- 
metric  system  in  (30)-(32). 

V.  SINGLE  STRIP  WITH  NORMAL  INCIDENCE 

When  a  plane  wave  has  normal  incidence  on  a  single  strip, 
the  scattered  field  can  be  expressed  as  follows: 


E‘ =  iEo  j  ^  cos  gye  </g 

■b 

(33) 

-£osgn(A.)  /■“  -,i»i  . 

Hy^  -  1  yj4  cosgye  '  dg 

(34) 

=g^-k^. 

(35) 

Let  us  express  the  cunent  density  on  the  strip  as 
series: 

a  Fourier 

N 

J  =  z  2  (.iifylb). 

0 

(36) 

(The  edge  mode  will  be  introduced  later.)  The  boundary  con¬ 
dition  on  tangential  //  leads  to 

4 

jf  yA  cos  gydg 

1  ^ 

j  -iO.Sjkfl/Eo)  2*/  cos  (Inylh),  —h<y<h 

1  0.  /i<l.v| 

(37) 

Multiplying  both  sides  of  (37)  by  cosgy  and  integrating  both 
.sides  over  the  range  0  <  y  <  «>,  we  obtain 

/EtJM’  ^ 

A  =  ~  *1^'  / 

ITT/:  0  0 

(38) 

r 

(/■/  =  ( 1  /  *'■)  /  cos  (/rr.t’/h)  cos  (gi  )  dy. 

(30) 
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The  boundary  condition  on  tangential  £'  leads  to 


1  A  cos  gy  dg 
^0 

=  -l. 

—h  <y  <  h. 

(40) 

Multiply  (40)  by  cos  (inylh)  and  integrate 
0<y  <h  to  obtain 

over  the  range 

/  AGidg== 
■b 

[-0.5, 

1  0. 

if  1  =  0 

if  1=^0. 

(41) 

From  (38)  and  (41), 

N 

2)  biZj,  =  K, 
1=0 

(42) 

Zit  =  (2/ki?w/ir) 

I  {Gfith)dg 
^0 

(43) 

where  K/  =  0  except  Fq  =  Eq. 

Now  let  us  express  the  current  density  on  the  strip  as  the 
sum  of  the  edge  mode  and  a  Fourier  series  as  in  (27).  Once 
again  we  obtain  (30)-(32)  with  Zn  given  by  (43). 


By  integrating  the  current  density  across  the  strip,  it  is 
found  the  current  is  /  =  —  h’Co  +  wc/v.  In  the  broadside 

direction  (0  =  0)  the  scattered  far-zone  field  is 

Es’  =  - (kvI/4)  . (44) 
ynkp 

For  normal  incidence,  the  backscatter  echo  width  is 

W,  =  2»rp  I  F*| V I  £o  1^  = - —  ,3-  • 

8fr|Eo  1 

VI.  NUMERICAL  RI.SULTS 

Fig.  2  illustrates  the  reflection  coefficient  versus  spacint; 
for  a  strip  grating  with  normal  incidence.  When  the  spacing  is 
less  than  the  wavelength,  these  data  show  satisfactory  agree¬ 
ment  with  the  measurements  of  Primich  1 10]  and  the  formu¬ 
las  of  Marcuvitz  (11). 

Fig.  3  shows  the  convergence  properties  of  our  calculations 
with  (22)  and  (30).  For  numerical  calculations  the  scries  in 
(32)  is  truncated  after  M  terms,  and  best  results  are  obtained 
if  M  is  greater  than  N  in  (30).  Good  results  are  obtained  with 
iW  =  20  and  =  5  for  0  <  w/X  <  0.5  or  /V  =  10  for  0.5  < 
w/X  <  I . 

Accurate  current  distributions  are  not  obtained  when  J  is 
expressed  as  a  Fourier  series  as  in  Section  III.  Suppose,  for 
example,  we  choose  N  =  20  in  (17)  and  (22).  Then  we  find 
the  calculated  current  distribution  has  20  prumtunced  rip¬ 
ples  across  the  strip.  These  false  ripples  are  not  observed  with 
the  edge-mode  formulation  in  Section  IV. 

For  a  strip  grating  with  normal  incidence  and  spacing  .< 
lesa  than  X,  it  is  found  (he  current  distribution  yt.rl  has  con¬ 
stant  phase.  Although  the  magnitude  of  J  varies  across  the 
strip,  the  phase  of  7  is  precisely  uniform.  This  is  observed 
for  all  values  of  the  strip  width  w.  Thus  yi.i  )  i|ualifirs  as  a 
“characteristic  mode"  (81,  (9)  of  the  grating,  since  this 


STRIP  SPACING  t/V 

Fig.  2.  Reflection  coefficient  of  perfectly  conducting  stiip  grating 
with  normal  incidence. 


t  ig.  3.  Cimvcrgcncc  curve  for  reflection  coefficient  of  strip  grating 
with  normal  incidence. 

current  distribution  generates  an  electric  fielil  which  also  has 
coiislani  phase  across  the  strip.  When  the  spacing  exceeds  X, 
the  phase  of  J  is  no  longer  constant. 

When  the  strip  width  is  just  slightly  smaller  than  the  spac¬ 
ing,  it  is  found  the  calculated  current  density  J  approaches  the 
physical  optics  approximation  7  =  2  //',  With  the  edge  mode 
formulation  of  Section  IV,  good  results  are  obtained  even 
with  large  spacing  such  as  r  =■  20X  on  =  50X.  When  the  spac¬ 
ing  is  large,  the  calculated  reflectiim  coefficient  agrees  closely 
with  Ihc  physical  optics  approximation  H  -  u  /.r. 

F'lg.  4  illustrates  the  backscatter  echo  width  of  a  single 
•strip  with  broadside  incidence.  Ihe  moment  method  (Section 
V)  shows  excellent  agreement  with  the  rigorous  solution 
Fig.  5  shows  Ihc  convergence  properties  of  out  calculations. 

Macrakis  (12)  presents  Ihe  rigorous  solution  foi  (he  strip 
ill  terms  of  Mathieu  functions.  We  used  Ihc  Mathicu  sub¬ 
routines  of  Hodge  1131  to  tabulate  the  rigorous  data  ,n 
I'ig.  4  and  Table  I.  Our  rigorous  data  agree  with  those  of 
Macrakis  in  Ihe  range  0..3I8  <  w/X  <  0  vs.s  where  tabulated 
Mathieu  fiinctiony  were  available  to  him.  However,  our  data- 
differ  greatly  from  his  (and  from  the  "exact"  data  ol  King 
and  Wu  I  14 1 1  in  the  region  u'/X  <  ').2.  I'or  narrow  strips 
with  ic/X  <  0.1 .  our  data  agree  with  calculations  based  on  the 
e<|uivalent  ruiiiid  wiie  with  radius  a  =  u'/4  (rather  than  a 
h72  as  used  in  1 1  2 1  ,ind  |  1 4 1 1. 

With  entire  basis  fiinclions.  convergence  greatly  improves 


15 


RICHMOND:  THEORY  OETM  SCATTERING 


887 


«• 


Fig.  4.  Broadside  backscatter  echo  width  of  perfectly  conducting 
strip. 
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Fig.  5.  Convergence  curve  for  broadside  backscatter  from  perfectly 
conducting  strip. 


when  the  edge  mode  is  included  among  the  basis  functions. 

It  remains  to  be  seen  whether  similar  improvement  can  be 

obtained  with  subsectional  bases  or  with  three-dimensional 

scattering  geometries  such  as  the  rectangular  plate. 
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On  the  Edge  Mode  in  the  Theory  of  Thick  Cylindrical 
Monopole  Antennas 

JACK  U.  RICHMOND.  FELLOW.  IEEE 


Abstract — The  Fourier  (ransform  is  employed  to  develop  Ihe  theory 
of  thick  cylindrical  monopole  antennas.  The  tubular  monupole  with  an 
open  end  is  considered  as  well  as  the  solid  monopole  with  a  flat  end. 
The  unknown  current  or  field  distribution  is  expanded  in  a  Fourier 
series,  and  Galerkin's  method  is  emptoyed  to  develop  simultaneous 
equations  for  the  Fourier  coefficients.  When  the  edae  mode  is  inciuded 
among  the  basis  functions,  it  is  found  this  greatiy  improves  the 
convergence  of  the  moment-method  calculations.  Numerical  data  are 
included,  showing  excellent  agreement  with  experimental  measure¬ 
ments  of  the  monopole  admittance. 

[.  INTRODUCTION 

Cinarsson  ( 1 1  presents  solutions  for  tubular  and  solid  cylin¬ 
drical-wire  antennas.  These  solutions  are  based  on  a  delta-gap 
model  of  the  generator  and  thus  arc  not  suitable  for  calcu¬ 
lating  the  susceptance.  This  objection  applies  also  to  the  for¬ 
mulations  of  Chang  12]  and  King  and  Wu  |31.  Otto  |4]  em¬ 
ploys  a  realistic  magnetic-frill  model  for  the  generator  and 
reduces  Einarsson's  formulation  to  one  equation  with  one 
unknown  (rather  than  eight  simultaneous  equations).  It  has 
not  been  demonstrated  that  Otto’s  approximations  are  appli¬ 
cable  to  electrically  thick  antennas,  however. 

Chang  (5)  presents  a  moment-method  solution  for  the 
thick  tubular  monopole  based  on  Halicn's  integral  equation 
with  an  extra  term  corresponding  to  a  radially  extended  source 
in  the  aperture  of  the  coaxial  feed  cable.  Thus  an  excellent 
solution  is  available  for  the  thick  tubular  antenna,  but  evi¬ 
dently  not  for  the  thick  solid  cylindrical  antenna. 

King  (111  presents  measured  admittances  of  thick  mono¬ 
poles,  theoretical  results  for  thick  tubular  monopoles,  and  a 
conection  term  for  higher  order  modes  in  the  coaxial  feed 
line. 

In  this  communication  we  present  moment-method  solu¬ 
tions  for  tubular  and  solid  cylindrical  antennas.  In  comparison 
with  previous  formulations,  ours  arc  straightforward.  It  is 
demonstrated  that  the  convergence  of  the  solution  is  greatly 
improved  when  the  edge  mode  is  included  among  the  basis 


Manuscript  received  January  J,  1980;  revised  April  29,  1980.  This 
work  was  supported  in  part  by  Contract  l)AA(;29-79<'-0082  between 
the  Army  Research  Otfrec  and  the  Ohio  State  University  Research 
I'oundation. 

The  author  is  with  the  I  lcctro.SVicnce  laboratory.  Department  of 
Electrical  I  ngineeting.  the  Ohio  Slate  University,  Columbus,  ()M  t  ,12 1 2. 
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functions.  Numerical  results  are  included,  showing  excellent 
agreement  between  the  calculated  and  measured  admittances. 

A  periodic  coUinear  array  of  cylindrical  wires  is  considered 
briefly.  When  the  incident  field  is  independent  of  the  coordi¬ 
nates  0  and  z,  an  interesting  property  is  observed.  I'or  all 
values  of  tlic  wire  radius  a,  the  wire  IcnKth  /.  and  llic. 
s  (with  z  <  X),  the  induced  currenl  dislrihulion  /(z)  has  pre¬ 
cisely  constant  phase  (i.e.,  the  phase  is  independent  of  z). 
The  electric  field  generated  by  this  current  also  has  constant 
phase  over  the  cylindrical  surface.  Thus  the  solution  (luatifies 
as  a  “characteristic  mode”  as  defined  by  Carbaez  (6|  and  Har¬ 
rington  and  Mautz  [7).  When  the  spacing  z  exceeds  X,  this 
constant-phase  property  is  lost. 


MACMETIC 


Fig.  I .  Center-fed  cylindrical  antenna. 

It  is  convenient  to  express  the  current  distribution  on  the 
antenna  in  a  Fourier  series  as  follows: 


II.  THE  CENTER-FED  TUBULAR  DIPOLE 

Fig.  1  illustrates  a  center-fed  tubular  dipole  antenna.  This 
dipole  is  equivalent  to  a  monopole  antenna  fed  via  a  coaxial 
cable  through  a  large  ground  plane.  The  tube  has  perfect  con¬ 
ductivity  and  infinitesimal  wall  thickness,  and  the  time  de¬ 
pendence  is  understood  for  the  electromagnetic  field.  The 
field  is  the  sum  of  the  incident  field  (£',  H‘)  from  the  mag¬ 
netic  frill  and  the  scattered  field  from  the  electric  surface  cur¬ 
rent  on  the  tube.  The  field  is  independent  of  the  angular  co¬ 
ordinate  0,  and  the  only  nonzero  field  components  are 
Ej,  and  The  scattered  field  may  be  expressed  as  follows; 

=  (Jklri)  f  AKo(ya)li(yfi)cos(g2)dg,  p<a  (I) 
■'o 

= -ijk/r})  f  AJoiya)Ki(yp)cos(gz)dg,  a<p  (2) 

A) 

T*=g*-k*  (3) 

where  k  =  wV/lie,  17  =  y/JiR,  In  *•’<1  £»  are  the  modified 
Bessel  functions,  .4  is  a  function  of  g,  and  (p,  0,  z)  denote  the 
coordinates  in  the  circular  cylindrical  system.  Except  for  the 
conducting  tube,  the  medium  is  homogeneous  with  parameters 
H  and  E.  The  scattered  field  component  E/  is  continuous 
across  the  boundary  at  p  =  a.  Setting  Ej{a,  z)  =  0  on  the  per¬ 
fectly  conducting  tube,  we  obtain  the  integral  equation  for 
the  tubular  dipole; 


j  y^lo(yo)flo(y«)  W*  Cf*)  dg  =  a). 


--/|<Z</|. 


(4) 

We  multiply  both  sides  of  (4)  by  cos  (inz/L  )  and  integrate 


over  the  region  0  <  z  <  X  to  obtain; 

f  y->*io(yay^o(y<>)G,dg  =  -v, 

(5) 

*0 

G,»  j 

r" 

f  cos  (gz)  cos  (ittz/f.)  dz 

(6) 

A 

> 

"'“i 

f  E,'(a,  z)  cos  (inz/L)  dz 

i 

(7) 

'0 


/(z)=  ^1>I  cos  (lifz/L)  (8) 

I 

where  /  is  a  positive  odd  integer.  The  scattered  field  must  sat¬ 
isfy  the  boundary  condition 


///(<» +.  r)-V(«-  0  =  -/z(0 
|/(zy(2irfl).  -h<z<h 
jo,  X  <  z  <  » 


•/z(a)  = 


(9) 

(10) 


where  denotes  the  sum  of  the  electric  surface-current  den¬ 
sities  J2  on  the  inner  and  outer  surfaces  of  the  tube.  From  ( I ), 
(2).  and  (9). 

(Ah)  cos  (gz)dg=JJ(z).  (ID 

We  multiply  both  sides  of  ( 1 1 )  by  cos  (g'z )  and  integrate  over 
the  region  0  <  z  <  “>  to  obtain 

A=—2jbfit.  (12) 

It  K  I 

In  (5)  let  us  replace  A  with  the  equivalent  quantity  from  ( 12). 
This  leads  to  the  following  system  of  simultaneous  linear  equa¬ 
tions; 

AT 

"^b,z„=v„  1=  1.3,5.  (13) 

r=i 

Z,I  =  —  I  y  lo(ya)Eo(ya)(Jfiidg.  (14) 

»r  *  '0 

The  impedance  matrix  ///  is  symmetric.  The  current  dis¬ 
tribution  represented  by  (8)  is  a  periodic  :vcn  function  with 
period  4X.  Since  each  term  vanishes  when  z  =  X,  (8)  repre- 
senls  a  current  distribution  which  vanishes  at  each  end  of  the 
tubular  dipole. 

Corresponding  to  the  TEM  mode  in  the  coaxial  feed  cable, 
let  the  frill  in  Fig.  I  have  the  following  magnetic  surface-cur¬ 
rent  density: 


k 

Jv 


r 


M, 


-»>^ti 

p  In  (b/a) 


(IS) 


where/.  =  2X. 


where  F|  g  denotes  the  generator  voltage  and  X  is  the  outer 
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radius  of  the  coaxial  cable.  For  p  <  a,  the  free  space  field  of 
the  magnetic  frill  is  (8) 

V  =  — I  (l/7j[A'o(7a)-Afo(7b)l 

vri  In  (b/a)  Jg 

'^l(7P)cos(gz)dj.  (16) 

From  (7)  and  (16),  the  excitation  voltages  for  the  center-fed 
tubular  dipole  are 

y 

*'i  ='  ,  /  l^o(7a)-Ao(7b)lfo(7‘i)G/dg.  (17) 

» In  (6/fl)  Jq 

The  admittance  of  the  tubular  dipole  is  determined  as  follows: 
y,,  =/(0)/^'ii  =(l/l'ii)2*l  (18) 

I 

Now  consider  the  scattering  problem  where  a  plane  wave 
has  broadside  incidence  bn  the  conducting  tube.  The  incident 
field  is 

E,‘=EQe~’'‘*  ^EoJo(kp).  (19) 

From  (7)  and  (19),  the  excitation  voltages  are 

(20) 

in 

where  n  =  (i  —  I)/2.  Of  course,  if  the  cylinder  has  a  large 
radius,  the  higher  order  modes  with  cos(m0)  behavior  cannot 
be  neglected. 

Finally,  let  us  consider  a  periodic  collincar  array  of  tubular 
dipoles.  As  in  Fig.  2,  let  s  denote  the  centcr-tu-ccnler  spacing. 
Equations  (7),  (8),  and  (13)  apply,  and 

•• 

Zt,  =  —  e„y„^Io(.y„a)Ko{ynn)Gi„G,„  (21) 

nks  „=o 

y„^  =  (2nnls)^  -  (22) 

Gin  —  I  cos(inz/L)cos(2nm/s)dz  (23) 

where  e„  =  2  except  Cq  =  1.  When  a  plane  wave  has  broadside 
incidence  on  the  array,  the  voltages  F,-  arc  given  by  (20). 
Equations  (13)  and  (20)-(23)  were  programmed  fur  the  digital 
computer.  When  the  spacing  s  is  less  than  or  cr/ual  to  A,  it  is 
found  the  complex  Fourier  coefficients  hi  all  have  the  same 
phase  angle.  It  follows  from  (S)  that  Ihe  pha.se  of  /(;)  is  inde¬ 
pendent  of  z. 

in.  EDC.K  MODE  FOR  THE  TUBULAR  DIPOLE 

Let  us  express  the  current  distribution  on  the  tubular 
dipole  as  the  sum  of  an  edge  mode,  a  feed  mode,  and  a  Fou¬ 
rier  series,  as  follows: 

Hz)  =  /,(r)  +  Ij{z)  +  ^c,cos(lnz/L).  ( 24 ) 

I 


I - - - ^ 


Mg.  2  Periodic  eolhneai  anas  ol’olindcrs 


The  idea  is  to  separate  out  the  rapidly  varying  components  of 
the  current  distribution.  Then  the  Fourier  series  may  converge 
rapidly  since  it  represents  a  smooth  function.  With  the 
moment  method,  one  must  invert  the  impedance  matrix  Z,j 
to  determine  the  Fourier  coefficients  and  the  antenna  admit¬ 
tance.  Thus  any  technique  that  reduces  the  size  of  the  matrix 
Zfi  implies  a  reduction  in  computational  expense  and  storage 
requirements  and  the  capabilily  of  solving  longer  dipole  an¬ 
tennas. 

In  the  vicinity  of  each  end  of  the  tubular  dipole,  it  is 
known  that  /(z)  must  behave  like  the  square  root  of  the  dis¬ 
tance  from  the  end  (9| .  Thus  a  suitable  edge  mode  is 

(z/h)^  =  CM  2d;  cos  ilnz/l)  (25) 

/ 

d,  =  (2/l)Jt0n/2)  (26) 

where  7;  denotes  the  Bessel  f'metion.  There  is  only  one  un¬ 
known  con.stanl  c\  in  the  edge  mode. 

The  current  function  Hz)  often  has  large  slope  in  the  vicin¬ 
ity  of  the  feedpoint,  which  in  the  present  case  is  at  the  center 
of  the  dipole.  This  component  of  Hz)' is  the  same  as  the  cur¬ 
rent  distribution  /plzl  on  an  infinitely  long  cylindrical  an¬ 
tenna.  Thus  a  suitable  feed  mode  is 

//(r)=  W  cos  (nz/L)  cos  (Inz/L).  (27) 

I 

I  he  facloi  cos  ()iz//.  I  is  included  so  the  feed  mode,  like  all 
lilt-  olliei  modes  in  ( 24).  will  vanish  at  the  ends  of  Ihe  dipole. 
The  current  distribution  on  Ihe  infinitely  long  cylinder  with 
magnelit  frill  excitation  is  |8)  : 

[”  1Ao(7b)  -Ao(76)]  cosCgOdg 

lo(z)  = - 1 - z - 

ijln(fe/fl)4  7^80(70) 


(28) 

An  efficient  subroutine  is  available  for  calculating  Ihe  function 
/o(r).  The  Fourier  coefficients  /;  for  the  feed  mode  are  calcu¬ 
lated  via  nunierical  integialioii.  T  he  feed  mode  contains  no 
unknown  constants. 

From  («).  (24).  (25),  and  (27), 

6;  =  c, +  e/vt/, +/,.  (29) 

In  (13)  we  replace  A;  with  the  equivalent  quantity  from  (29) 
to  obtain  the  following  system  of  simultaneous  linear  equa¬ 
tions: 

A 

'^c,z,;  =  y;,  1  =  1,3.  ■•,  A'  (30) 

/=! 
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y/ 

/=i 

1<N 

Ziin'  =  2  dl^il 
1=1 


(31) 

(32) 

(33) 


//«  =  (jkh)  2 


Bil  I  (I'/P) 


1=0  tMyfi) 


-  (/fc/Ti)  S  ~~ cos  Onnzls),  p>a 
n  =  0  ^o(‘l'n“) 


2Cn 

~  ^1^ III  ^o(2f|0)(ln 

7nT  / 

f‘‘ 

Fi„=  I  cos(lm'/(i)cos(2iim/s)iJ!’. 


The  fielil  2)  must  match  the  field  r)  over 

the  aperture.  I'tiis  leads  to 

— — ; - =  X  (7,, <!)/•/„ 

n  Ao(7h«)  T 

«<7t/o(7/«) 

From  (37)  and  (39)  we  ohlain  the  simultaneous  linear  equa¬ 


tions: 

N 


'^B,Zi,=  Vi,  1  =  0.  1.-.  A' 


/«0 


7n''^o(7na) 


where  I'j  and  Z,7  are  given  by  ( 14)  and  (17).  A.s  Af  is  increased, 
the  admittance  converges  much  mure  rapidly  with  (30)  than 
with  (13).  The  impedance  matrix  Z,/'  is  not  symmetric.  Al¬ 
though  most  of  the  constants  c/  are  Fourier  coefficients,  the 
last  one  c/v'  represents  the  amplitude  of  the  edge  mode.  The 
quantities  d/  and  //  are  known  constants. 

In  the  plane  wave  scattering  problem,  it  is  better  to  delete 
the  feed  mode.  In  this  case,  K,  and  F/  are  given  by  (20). 

IV.  PERIODIC  COLLINEAR  ARRAY  OF  SOLID  DIPOLES 

Let  us  consider  a  periodic  collinear  array  of  solid  perfectly 
conducting  cylindrical  dipoles  with  flat  ends,  as  shown  in  Fig. 
2.  The  incident  field  (£',  //')  is  considered  to  be  periodic  also 
with  period  s,  and  the  source  is  located  in  the  exterior  region 
where  p  >a.  The  field  may  be  expressed  as  follows: 


yW/i(7/a)6„ 

Zit  = - - - -  +(2/*/r)2v 

errMyfi)  „ 


i{y  In 

ynFoiyno) 


(40) 


(41) 


(4:> 


cos  (Jnz'ld),  p<a  (34) 


(35) 


=  (jklr\)'^  Ent\(ynP)f^o%i2m[zl$\  p<.a  (36) 

n 

where  r'  =  2  —  A.  It  is  easy  to  verify  that  (34)  represents  a 
Held  with  Kp  vanishing  over  the  flat  ends  of  the  cylinders  at 
a'  =  0  and  z  =  2d.  Equation  (34)  is  valid  only  in  region  / 
where  h  <  z  <  h  +  2d.  Of  course,  the  field  vanishes  in  (he  in¬ 
terior  of  each  solid  cylinder. 

The  field  E,(a+,  z)  must  vanish  on  the  surface  of  the  con¬ 
ducting  cylinders,  and  it  must  match  Ei{a—,  z)  over  the  aper¬ 
ture  between  cylinders.  Thus 


As  usual,  we  assume  the  incident  field  (and  the  Fourier  coeffi¬ 
cients  £„)  is  known.  The  B/  are  determined  from  (40)  by 
matrix  inversion  or  equivalent  techniques.  Then  the  A„  are 
obtained  from  (37).  Although  the  details  cannot  be  ptcsenled 
here,  one  can  then  determine  the  current  distribution  and  the 
admittance  of  the  solid  cylindrical  antenna.  The  aperture  field 
can  be  expressed  as  the  sum  of  an  edge  mode  plus  a  well- 
beha>'ed  residue  with  techniques  analogous  to  those  presented 
for  the  tubular  cylinder. 

V.  NUMERICAL  RESULTS 

Let  us  discuss  first  the  numerical  results  for  the  tubular  an¬ 
tenna.  Fig.  3  illustrates  the  convergence  curves  for  the  admit¬ 
tance  of  a  tubular  dipole.  M  denotes  the  size  of  the  matrix  Z,( 
(i.c.,  the  number  of  simultaneous  linear  equations)  In  the 
moment-method  calculations.  A  dramatic  improvement  is  to 
be  noted  in  Fig.  3  when  the  current  dislrihuliou  is  represented 
as  the  sum  of  the  edge  mode,  the  feed  mode,  and  a  Fourier 
series.  Convergence  is  obtained  with  M  =  6,  whereas  M  must 
exceed  40  when  the  current  is  represented  by  a  Fourier  series 
alone. 

Fig.  4  illustrates  the  admittance  versus  length  for  a  tubular 
monopole.  Our  calculations  show  excellent  agreement  with  the 
experimental  measurements  of  Holly  110|.  In  these  calcula¬ 
tions  accuracy  was  more  important  than  computational  econ¬ 
omy.  Therefore,  we  used  Af  =  10  for  Fig.  4.  In  the  infinite 
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MATRIX  SIZE  M 

Fig.  S.  Convergence  curve  for  admiltance  of  solid  center-fed  cylin¬ 
drical  antenna  with  flat  ends. 


series  in  (31)  anil  (33),  let  LL  denote  (he  number  of  turnis 
actually  summed.  We  used  l.l,  =  20  for  Fig.  4. 

Fur  (he  solid  cylindrical  antenna.  Fig.  5  sliows  substantial 
improvement  in  convergence  when  the  edge  mode  is  included 
among  the  ba.si.s  functions.  Fig.  6  illustrates  the  admittance 
versus  spacing  for  a  periodic  collinear  array  of  solid  dipoles. 
Usually,  one  is  interested  in  a  single  dipole  rather  than  an 
array.  The  admittance  converges  rapidly  with  K  if  s/X  =  K  I 
1/2  where  K  denotes  an  integer.  Thus  the  array  admittance  is 
essentially  equal  to  the  single-dipole  admitlance  when  s/X  = 
3.S  in  Fig.  6. 

Fig.  7  illu.strarcs  the  adrnittjncc  vcr5ti5  length  for  j  solid 
cylindrical  monopolc.  Again,  the  ealculutions  show  excellent 
agreement  with  experimental  measurements.  In  these  caleula- 
lions  we  used  Af  -  14  and  !J,  =  10,  In  the  ealcutations  for  Fig. 
7  wc  lei  i/X  =  5.5, 


I'ig.  6.  Admittance  versus  spacing  for  periodic  collinear  array  of  solid 
center>fed  cylindrical  antennas  with  flat  ends. 
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evaluate  Ihe  eurrenl  are  si^en.  The  transient  current  induced 
on  the  wire  is  compared  to  the  current  computed  by  assuming 
quasi-steady-slale  theory. 

The  interpretations  and  restrictions  of  all  the  given  solu¬ 
tions  are  discussed 
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Simple  rigorous  expressions  are  developed  for  the  Helds  of  an  electric  source  current  distribution 
that  is  separable  in  cartesian  coordinates.  The  medium  is  unbounded  free  space,  and  the  time-hartnonic 
source  current  exists  in  the  slab  region  s^<  z  <  i,. 


INTRODUCTION 

Many  textbooks  present  the  theory  of  separable 
fields  in  a  source-free  homogeneous  region. 
However,  they  do  not  relate  these  fields  to  the 

A(x,  y,  z) 

(6) 

sources  or  determine  the  fields  in  the  source  region. 

If  X  and  Y  satisfy  (2)  and  (3),  it  is  shown  in  Appendix 

These  relations  and  expressions  are  useful,  for 

1  that 

example,  in  the  study  of  propagation  in  periodic 
media,  radiation  from  a  periodic  array  of  sources. 

f  f 

and  scattering  from  a  periodic  array  of  metallic 

1  X{x')Y(y')  —  dx'dy' 

or  dielectric  elements. 

J  J 

•ea 

A 

THEORY 

0) 

Consider  a  time-harmonic  electric  current  density 

7?  +  y'j  +  y!  =  =  -o>V* 

(8) 

with  the  form 

From  (6)-(7),  the  vector  potential  is 

J  =  (£A  +yB  +  zC)  X(x)  Yiy) F(z)  (I) 

A  =  +  yfl  +  zC)XYZI(lyy) 

(9) 

where  A,  B,  and  C  are  arbitrary  constants,  the 
time  dependence  e^**'  is  suppressed,  F(z)  is  an 

1 

F(r')e’”'"‘ 'dr'  = 

(10) 

arbitrary  function,  and 

0  r  <  ri 

• 

X'^y^X  (2) 

y-yXY  (3) 

P  s  * 

1  F{z' )*’’'•' dz'  r,  <z<r, 

J  «i 

(H) 

Primes  denote  derivatives,  and  y,  and  -y,  arbi¬ 

f  *1 

trary  complex  constants.  Suitable  functions  are 

\  Fiz'ye-"-  dz'  z,<i 

=  cosh  7,x,  sinhyiX,  *■*'',  e"''*  (4) 

0  *1  <  * 

Y{y)  =  cosh y^y,  sinh y,>,  e''’’,e~'^’'  (5) 

j  F(z')«-''’  dz’  z,<z<z. 

or  linear  combinations  of  these  functions. 

Q‘  • 

(12) 

Let  the  source  J  exist  only  in  the  slab  region 

f  ** 

r,  <  r  <  r,.  If  the  medium  everywhere  is  free 

\  dz*  /</, 

space,  the  vector  potential  is  given  as  follows: 

J  '1 

Z"-yiZ-2y,f 

(13) 

Ciwulilll  e  INI  ky  dM  XwwtcM  OMfkyUnl  Uaim. 

The  magnetic  field  intensity  is  found  as  follows; 

Ptpur  numbur  IS09I4.  1299 

004S-64(M/8I  / 1 1 12-09I4$01 .00 
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II  -  7  »  A  n4j 

H^^^iCXY'Z-  BXYX‘)I(2'1,)  (15) 

H,  =  (AXYZ'  -  ex'  rZ)/(2-y,)  (16) 

H.=^(BX' YZ- AXY‘Z)l(2-i^)  (17) 

The  electric  field  intensity  is  given  by 

E  =  (VxH-J)/(7W)  (18) 


»  l(fc'  +  'i\)AXYZ  +  BX'  i  'Z  +  ex'  YZ'  \  D 

(19) 

E,  =  [AX'  Y'Z+(k^  +  yl)BXYZ+eXY'Z']D 

(20) 

£.=  [AX'  YZ'  +  BXY'Z' 

+  (k^ +  y\)eXYZ-2y^eXYF\D  (21) 

D=l/(27W-y,)  (22) 

In  (21)  the  last  term  vanishes  if  z  <  z,  or  2  > 
With  the  aid  of  (1H3),  (8).  (9),  and  (13).  it  is 
easy  to  show  that  our  vector  potential  A  satisfies 
the  appropriate  wave  equation 

V’A  +  fc'A=-J  (23) 

Since  our  solution  in  (15)-(17)  and  (19)-(2I)  is 

derived  via  (14)  and  (18),  it  is  not  difficult  to  show 
that  is  satisfies  all  four  of  Maxwell’s  differential 
equations.  In  Appendix  2  it  is  shown  that  it  also 
satisfies  the  required  boundary  conditions  at  z  = 
z, .  The  radiation  conditions  at  z->  oo  will  be  satisfied 
if  we  choose  the  proper  root  for  -y,  . 

Theorems  are  available  that  make  it  possible  to 
use  the  above  free-space  fields  even  when  the 
medium  is  not  free  space.  For  example,  suppose 
one  has  a  periodic  array  of  perfectly  conducting 
wires  surrounded  by  free  space.  It  is  well  known 
that  the  field  will  not  be  disturbed  when  the  wires 
are  removed  if  the  surface-current  density  is  main¬ 
tained  in  free  space.  Next,  suppose  one  has  a 
periodic  array  of  dielectric  rods  surrounded  by  free 
space.  Without  disturbing  the  field,  one  may  remove 
the  rods  while  maintaining  the  polarization-current 
density  J  =2  y'oife,  -  e)E  in  free  space.  If  the 
permeability  p,  of  the  rods  differs  from  ix,  then 
the  magnetic  current  density  M  =  yuifiJi,  -  |i)ll 
is  required  also. 

Suppose  one  has  a  periodic  array  of  wires  or 
rods  in  the  region  z,  <  z  <  z,  and  the  surrounding 
medium  is  free  space.  Then  one  may  write  a  suitable 


expression  for  the  current  density  J(x,  y,  z)  and 
apply  the  Fourier  transform  with  respect  to  the 
X  and  y  variables.  The  resulting  expression  for  J 
will  be  a  summation  of  separable  current  functions. 
Each  term  will  have  the  form  of  (1)  with  a  subscript 
n  inserted  on  the  quantities  A,  B,  C,  X,  and  Y. 
Then  the  field  of  the  periodic  array  is  given  by 
(IS)-(2I)  with  insertion  of  subscripts  and  summa¬ 
tions. 

CONCLUSION 

Simple  rigorous  expressions  are  developed  for 
the  fields  of  an  electric  source  current  distribution 
that  is  separable  in  cartesian  coordinates.  The 
medium  is  unbounded  free  space,  and  the  time-har¬ 
monic  source  current  exists  in  the  slab  region  z, 
<  Z  <  Zj. 

APPENDIX  1: 

EVALUATION  OF  THE  DOUBLE  INTEGRAL 

The  double  integral  in  (7)  can  be  evaluated  most 
readily  with  a  technique  presented  by  Harrington 
[1961].  In  this  method,  one  solves  a  radiation 
problem  by  two  different  methods  and  equates  the 
two  solutions.  To  apply  this  method,  let  us  consider 
an  infinite  sheet  of  electric  surface-current  density 
J,  =  xX(x)Y(y)  located  on  the  plane  surface  z 
—  z'  and  radiating  in  free  space,  where  X  and  Y 
satisfy  (2)  and  (3).  From  the  boundary  conditions 
and  symmetry  considerations,  H^  =  ±XY / 2  when 
z  — »  z' .  It  is  easy  to  deduce  that 

//,=  -(I/2)2r(x)y(.v)sgn(z-z')e (24) 
and  therefore  the  vector  potential  is 

A.  =  2f(x)r(>-)e '  7(27,)  (25) 

But  the  vector  potential  is  also  given  by 

-<.=  [(  X(x’)Y(y’)'—<ix'dy'  (26) 


By  equating  these  two  expressions  for  A , ,  we  obtain 

(7). 

An  alternative  derivation  of  (7)  is  available  if 
7i  =  if  I  and  7j  =  ygj  in  (2)  and  (3),  where  /,  and 
g2  are  real  constants.  The  following  expansion 
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t. 


(2«) 

is  developed  by  Tyros  [1969].  We  substitute  (27) 
into  (6),  interchange  the  order  of  the  integrations, 
and  note  that 

j  dx'  =  27r8(/  +  /,)  (29) 

J  e^'»**>'''>'  =  2ir8(«  +  g,)  (30) 

This  procedure  leads  again  to  (7)  and  (8). 


APPENDIX  2: 

THE  BOUNDARY  CONDITIONS 


In  order  to  be  correct,  the  solution  presented 
in(13)-(21)must  satisfy  certain  boundary  conditions 
at  the  surfaces  z  =  z,  and  z  =  Zj.  For  the  sake 
of  brevity,  let  us  consider  only  the  boundary  at 
z  Zj.  From  (IS),  (16),  (19),  and  (20),  the  tangential 
field  components  and  are  continuous 

across  this  boundary  if  Z  and  Z'  are  continuous. 
We  recognize  these  as  the  appropriate  boundary 
conditions  for  a  surface  that  does  not  contain  a 
sheet  of  surface  current  density.  From  (10)-(12), 


Z(z). 


+  *' 


Z'(z)  -  — yj« 


(  F(z')e’'‘' dz' 

J«l 

j  f(z' )«-’•'  dz' 

*(  F{z')e-^‘  dz' 


(31) 


F(z’)e-^‘ 

'  dz' 

(32) 

Zj,  and 

r)  -  «-•*  J 

•j 

f(z')«’>'  dz’ 

•l 

f  •> 

(33) 

1  F(z')e"’  dz' 

(34) 

for  z  >  Zj.  If  the  function  F(z)  is  well  behaved 
at  z  -  Z],  a  study  of  (3l)-(34)  shows  that  Z  and 
Z'  are  continuous  across  the  boundary  at  z  >  z,. 
Thus  the  tangential  field  components  satisfy  the 
required  boundary  conditions.  From  (17)  it  also 
follows  that  H,  is  continuous  as  required  in  the 
given  circumstances. 

On  the  surface  z^  Zj,  the  surface  charge  density 

is 
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p.=  CXYF(z,)/(jio>  (35) 

The  required  boundary  condition  on  the  normal 

component  of  E  is 

p.  =  tE‘.(z,)  -  tE'Xz,)  (36) 

where  z  >  z,  in  region  I  and  z,  <  z  <  Zj  in  region 
n.  In  (21)  for  £,,  all  the  terms  arc  continuous 
except  the  last  one.  Since  £(z)  vanishes  in  region 
1,  it  is  evident  that  our  solution  does  satisfy  (36). 
Finally,  let  £(z)  consist  of  a  unit  impulse  function 


£(z)  =  S(z  -  z").  From  (10)-(12), 

2(2)  =  (37) 

Z'  (2)  -  -7,  sgn  (2  -  2*)e-’'>''-'‘'  (38) 

Thus  Z  is  continuous  across  the  boundary  at  z  = 
z".  From  (15)-(22), 

H\  -  H“  =  BXY  (39) 

H',- H],' = -AXY  (40) 

E[  -  E"  -  -CX’  Y/{j,oE)  (41) 

El  -  E';  = -CXr /U<»t)  (42) 

where  z  >  z*  in  region  I  and  z  <  z*  in  region 
II.  The  electric  surface  current  density  at  z  =  z* 
is 

J.  =  (id +>g)2fy  (43) 


From  the  equivalence  principle  of  Mayes  [I9S8], 
the  electric  current  density  can  be  replaced  with 
an  equivalent  magnetic  surface  current  density  given 
by 

M.  -  (-xCXY’  +  yCX'  Y)l(j^c)  '  (44) 

It  is  easy  to  show  that  (39)-(42)  satisfy  the  required 
boundary  conditions  at  the  surface  z  =  z",  which 
contains  electric  and  magnetic  surface  current 
densities. 
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Many  buried  scattering  objects  of  interest  lake  the  form  of  two-dimensional  geometries.  This  in¬ 
cludes  scattering  from  utility  lines,  tunnels,  and  geological  structures  such  as  fault  lines.  Radar  systems 
used  to  detect  these  objects  commonly  use  antennas  that  are  at  least  comparable  in  sire  to  the  dep,h  of 
the  target.  Thus  the  target  and  the  antenna  do  not  satisfy  far  zone  conditions,  and  the  radar  range 
equation  is  not  applicable.  Consequently,  the  usual  separate  analyses  of  range,  antenna,  and  target  are 
not  possible.  In  a  previous  paper  a  method  was  outlined  that  permitted  the  antenna  properties  and  the 
scattering  properties  of  a  two-dimensional  target  to  be  treated  separately  for  the  case  of  a  linear  electric 
or  magnetic  dipole  source  parallel  to  the  axis  of  the  two-dimensional  scatterer.  This  involved  com¬ 
puting  the  received  voltage  for  an  antenna  located  at  the  image  position,  i.e.,  at  twice  the  target  range, 
and  computing  the  backscattered  fields  for  an  electric  or  magnetic  line  source  at  the  position  of  the 
transmitting  antenna.  This  model  has  also  been  applied  approximately  to  a  video  pulse  radar  with  an 
orthogonal  dipole  antenna  system.  It  would  also  be  applicable  to  large  loop  antennas  quite  commonly 
used  in  geophysical  explorations,  as  will  be  discussed.  The  primary  goal  of  this  paper  is  to  discuss 
additional  scattering  analyses  that  could  be  used  to  extend  the  previous  results.  The  major  thrust  then 
is  to  generate  solutions  for  the  scattering  attenuation  function  (SAP),  which  has  the  form  £*/£'  or 
H‘/H'  where  the  £'  (or  H‘)  are  the  electric  (or  magnetic)  fields  of  an  electric  (or  magnetic)  line  source  at 
the  image  position  and  the  £’  (or  H")  are  the  respective  scattered  fields.  Eigenfunction  solutions  have 
been  used  to  obtain  the  SAP  for  circular  cylindrical  geometries  to  represent  pipes  and  tunnels.  Moment 
method  solutions  have  been  applied  to  perfectly  conducting  wires  with  and  without  a  dielectric  sheath. 
Moment  method  solutions  have  been  applied  to  noncircular  penetrable  bodies  using  the  polarization 
currents  to  represent  the  unknowns.  Such  solutions  were  developed  for  a  line  source  above  dikes  by 
Parry  and  Ward  (1971)  and  Hohman  (1975)  in  the  early  seventies.  Such  solutions  can  be  made  for 
frequencies  that  include  several  target  resonances  for  potential  target  identification.  These  solutions  can 
possibly  be  extended  to  include  fault  lines,  joints,  etc.,  provided  their  electrical  properties  can  be 
estimated  by  using  some  of  the  concepts  involved  in  the  hybrid  geometrical  theory  of  diffraction- 
moment  approach.  The  methods  of  the  modified  geometrical  optics  could  also  be  applied  to  obtain 
scattered  fields  at  higher  frequencies.  These  and  other  potential  approaches  will  be  discussed. 


INTRODUCTION 

Scattering  from  buried  objects  differs  from  that  in 
free  space  in  several  significant  respects.  First,  there  is 
automatically  a  two-layer  geometry  since  an  earth- 
free  space  interface  is  present  with  the  probes  being 
located  at  the  interface  or  in  free  space.  This  can  be 
accounted  for  when  the  radiators  are  infinite  line 
sources  or  infinitesimal  electric  or  magnetic  (small 
loops)  dipoles.  However,  the  analysis  becomes  more 
complex  for  finite  antenna  geometries  and  generally 
involves  substantial  computation  time.  In  this  paper 
we  consider  only  the  scattering  from  lossy  objects 
immersed  in  lossy  homogeneous  media. 
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Actually,  P  irry  and  Ward  [1971]  partially  remove 
the  half  space  in  part  of  their  analyses  of  the  scatter¬ 
ing  from  dikes.  They  included  the  half  space  in  the 
field  incident  on  the  dike  but  not  in  the  fields  scat¬ 
tered  from  the  dike.  Comparison  with  the  solution 
when  the  half  space  is  completely  included  gave  vari¬ 
ations  of  the  order  of  10%.  This,  of  course,  is  trivial 
in  comparison  to  variations  in  geometry  for  ‘real 
world'  situations. 

Second,  losses  in  the  ambient  medium  occur  as  the 
wave  propagates  through  the  medium  and  over  the 
target. 

Third,  many  targets  which  are  to  be  detected  using 
a  three-dimensional  antenna  system  are  in  essence 
two-dimensional,  e.g.,  utility  lines,  tunnels,  faults,  and 
joints.  Consequently,  it  becomes  impractical  to  de¬ 
scribe  the  target  in  terms  of  far  field  concepts  such  as 
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Fig.  1.  Introducing  the  scattering  attenuation  function  of  a  cylin¬ 
der  into  the  radar  model. 


echo  area  (or  radar  cross  section)  or  echo  width.  One 
definition  that  partially  alleviates  this  difllcult  is  that 
of  the  scattering  attenuation  function  (SAF)  IBurrell 
et  ai,  1979].  This  parameter  will  be  defined  in  the 
next  section.  H  reduces  the  scattering  problem  to  that 
of  a  line  source  in  the  vicinity  of  the  two-dimensional 
target. 

Ward  [1980]  has  summarized  the  significant  devel¬ 
opments  of  electromagnetic  scattering  analysis  as  it 
pertains  to  the  lossy  earth.  In  his  notation,  the  ap¬ 
proach  used  herein  reduces  the  three-dimensional 
source  two-dimensional  inhomogeneity  to  a  two- 
dimensional  source  two-dimensional  inhomogeneity. 

Stoyer  and  Greenfield  [1976]  have  previously 
considered  the  three-dimensional  source  two- 
dimensional  target  geometry  in  the  presence  of  an 
interface.  However,  they  used  a  point  (infinitesimal) 
source  and  Fourier  transform  techniques.  The  ap¬ 
proach  to  be  introduced  here  makes  it  possible  to  use 
a  finite  length  source. 

Pridmore  [1978]  also  worked  the  three- 
dimensional  source  three-dimensional  target  con¬ 
figuration.  Richmond  [1978]  has  treated  this  problem 
where  the  three-dimensional  source  was  an  infinitesi¬ 
mal  magnetic  dipole  and  the  three-dimensional 
target  was  a  conducting  sphere  as  modeled  by  a  wire 
grid.  Actually,  the  wire  configuration  could  be  ar¬ 
ranged  to  fit  any  conducting  shape.  It  is  observed 
that  most  analytic  results  obtained  to  date  in  the 
geophysical  literature  are  for  scatterers  in  the  so- 
called  Rayleigh  region,  i.e.,  below  resonance,  whereas 
the  results  discussed  herein  include  the  resonance 
region. 

Most  of  the  work  at  the  ElectroScience  Labora¬ 
tory  has  been  directed  toward  physically  small  tar¬ 


gets  nearer  the  surface  at  higher  frequencies.  The 
medium  in  this  case  is  a  lossy  dielectric  in  contrast  to 
the  conducting  medium  used  throughout  the  geo¬ 
physical  literature.  Further,  many  of  the  targets  of 
interest  are  hollow  cylinders,  e.g.,  plastic  pipes  and 
tunnels.  Thus  one  must  be  cautious  in  extending  all 
of  these  results  to  the  low-frequency  conducting 
media.  Nevertheless,  if  one  uses  such  parameters  as 
index  of  refraction  and  impedance,  then  this  difficulty 
should  be  averted.  There  is  one  other  significant  ad¬ 
vantage  to  working  with  shallower  targets;  that  is  the 
ability  to  be  able  to  measure  the  scattered  fields  and 
to  compare  them  with  computed  values  \_Davis, 
1979]. 

The  goal  of  this  paper  is  to  review  the  results 
achieved  thus  far  for  evaluating  the  SAF  and  to  dis¬ 
cuss  techniques  for  extending  this  effort.  To  this  end, 
we  first  discuss  the  means  that  have  been  used  of 
evaluating  scattered  fields  from  penetrable  two- 
dimensional  scatterers  in  a  lossy  medium.  We  shall 
then  review  briefly  the  pertinent  solutions  for  the 
lossy  media  and  finally  indicate  how  they  might  be 
extended  to  the  geometries  of  current  interest. 

The  next  sections  of  this  paper  focus  attention  on 
the  analyses  we  have  used  to  evaluate  SAF's  for  lossy 
media.  Later  sections  then  introduce  techniques  that 
have  been  used  for  dill'erent  problems  that  could  be 
adapted  to  this  purpose. 

THK  SCATIl-RtNO  ATri;Nt)ATtON  l  UNt  TION 

The  scattering  attenuation  function  is  defined  as 

SAF  =  £•/£■  or  H’/H‘  or  V/U'  (I) 

where  U’  is  the  scattered  field  from  the  two- 
dimensional  target  at  the  line  source  position,  U'  is 
the  radiated  field  of  the  line  source  at  the  image  posi¬ 
tion  as  defined  in  Figure  I,  and  E  and  H  represent 
electric  and  magnetic  field  intensities  respectively. 

The  field  parameter  U  can  represent  either  the 
electric  or  the  magnetic  field  intensity  depending  on 
whether  the  source  is  an  electric  or  a  magnetic  line 
source,  it  is  also  possible  to  use  line  dipole  sources. 
Hohnum  [1971]  introduced  a  normalization  that  in¬ 
volved  the  magnetic  fields  in  free  space,  i.e..  Ampere's 
law.  The  normalization  introduced  here  tends  to  cast 
the  electrical  parameters  in  a  more  usable  form  as  we 
shall  soon  sec. 

The  basic  model  shown  in  Figure  1  introduces  a 
three-dimensional  antenna  at  the  image  position.  If 
the  source  antenna  is  also  the  receiving  antenna,  then 
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the  normalized  received  voltage  is  approximately 

=  (^V^V)SAF  (2) 

The  ray  spreading  in  the  plane  of  the  paper  of  Figure 
I  is  contained  in  the  factor  .  This  also  contains  a 
spreading  factor  in  the  perpendicular  plane  which  is 
canceled  by  the  ray  spreading  contained  in  the  factor 
£'  and  is  then  replaced  by  the  spreading  factor  in¬ 
herent  in  E^.  If  orthogonal  transmit-receive  dipole 
antennas  are  used  rotated  by  45°,  then 

F,  =  i(F;/F^)SAF  (3) 

Strictly,  (3)  is  valid  only  when  the  scattered  field  has 
a  preferred  polarization  independent  of  the  polariza¬ 
tion  of  the  incident  field.  This  can  be  modified  to 
account  for  the  case  where  the  scattered  fields  are 
simply  different  for  different  incident  polarizations 
[Burrell  et  al.,  1979],  This  model  has  been  confirmed 
by  comparison  with  a  moment  method  solution 
using  a  1000-m-long  conducting  cylinder  as  a  target. 
The  results  for  a  video  pulse  radar  agree  within  0.6 
dB  where  the  received  peak  signals  were  143  dB 
below  the  transmitted  signal  and  the  waveforms  are 
identical.  The  reader  is  referred  to  Burrell  et  al. 
[1979]. 

The  ratio  V^/yr  can  now  be  computed  separately, 
and  there  is  no  need  to  restrict  the  actual  antenna  to 
either  an  infinitesimal  dipole  or  a  line  source. 

it  is  apparent  that  the  SAP’s  are  a  function  of 
depth  d  and  the  electrical  parameters  of  the  medium. 
However,  these  effects  can  be  easily  estimated  with 
reasonable  accuracy  without  repeating  the  entire 
analysis.  Thus  the  SAF  of  a  specific  scatterer  can  be 
estimated  reasonably  accurately  (within  several 
decibels)  as  the  range  and  electrical  properties  of  the 
media  vary. 

The  scattering  attenuation  function  should  not  be 
seriously  changed  if  a  ground-air  interface  is  included 
in  the  model,  since  it  would  influence  the  scattered 
field  and  the  image  field  in  the  same  manner.  Two 
conditions  must  be  met  in  order  for  the  above  as¬ 
sumption  to  be  valid.  First,  there  is  no  significant 
multiple  target-interface  interaction,  and  second,  the 
lateral  extent  of  the  cylinder  target  is  sufficiently 
small  so  that  lateral  waves  on  the  surface  arc  not 
significant.  This  could  be  relaxed  to  the  extent  that 
the  extreme  rays  arc  not  incident  at  the  surface 
beyond  the  critical  angle. 

The  guidelines  to  follow  have  been  developed  for 
the  case  where  the  pertinent  scattering  mechanisms 
involve  internal  reflections.  The  scattered  fields  will 


then  be  proportional  to  the  Fresnel  reflection  coef¬ 
ficient. 

Guideline  I.  The  level  of  the  SAF  in  the  target 
resonance  region  is  changed,  as  the  electrical  pa¬ 
rameters  of  the  ambient  medium  are  changed,  by  the 
change  in  the  Fresnel  reflection  coefficient. 

Guideline  2.  The  changes  in  the  electrical  proper¬ 
ties  of  the  ambient  medium  alter  the  loss  contained 
in  the  SAF  by  an  amount  proportional  to  twice  the 
distance  from  the  leading  portion  of  the  (penetrable) 
target  to  the  reference  plane.  This  is  because  the  ref¬ 
erence  field  U‘  is  modified  by  this  change. 

Guideline  3.  As  the  depth  d  is  changed,  the  SAF  is 
changed  according  to  the  approximate  relation 


where  R  is  the  radius  of  curvature  at  the  first  station¬ 
ary  point  on  the  target.  This  is  based  on  geometrical 
optics  and  assumes  cither  that  this  specular  scatterer 
is  dominant,  which  is  usually  the  case  for  a  lossy 
medium,  or  that  R  «  d. 

These  guidelines  have  been  used  for  circular  tun¬ 
nels  [Burrell  et  al.,  1979]  and  have  worked  very  well 
for  the  relations  of  the  SAF  as  a  function  of  media 
and  depth.  Thus  the  SAF  parameter  is  somewhat 
more  general  than  it  would  originally  have  appeared 
and  in  practice  is  a  more  useful  concept  than  the 
conventional  echo  width. 

Finally,  one  observes  that  the  interface  could  have 
been  incorporated  in  the  model  of  Figure  1.  There 
would  be  no  substantial  difficulty  in  analyzing  the 
geometries  of  Figures  In  and  Ih.  However,  if  the  an¬ 
tennas  of  Figure  Ic  arc  of  finite  extent,  then  very 
large  computer  running  times  are  encountered  [Bur¬ 
rell  and  Peters,  1979].  More  recently,  using  tech¬ 
niques  developed  for  striplinc  antennas,  this  difficulty 
has  been  largely  eliminated  [Ucunoy/u  et  al.,  1979], 
but  this  has  not  been  incorporated  in  these  results. 

MODAL  SOLUTIONS 

The  sphere  and  the  infinite  cylinder  represent 
clii.ssic  shapes  in  the  study  of  electromagnetic  scatter¬ 
ing.  Prior  to  the  advent  of  large  digital  computers  it 
was  most  difficult  to  evaluate  numerically  the  scat¬ 
tered  fields  of  such  lossy  targets  even  when  the  ambi¬ 
ent  medium  was  free  space.  It  is  not  surprising  that, 
as  Parry  and  Ward  [1971]  observe.  D'Yakonov 
[1959(1,  6]  ‘has  published  a  solution  to  scattering 
from  a  circular  cylinder  and  a  sphere  in  a  conducting 
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Fig.  2.  Comparison  of  scatlering  altenuaiion  function  fc'n  for 
an  electric  line  source  for  t-m-radiu.s  tunnel  for  different  depths  d 
and  conductivity  a  shows  that  the  primary  effect  of  these  pa¬ 
rameters  is  to  change  the  amplitude  of  the  curve. 

half-space,  and  we  have  yet  to  witness  numerical  re¬ 
sults  for  either  of  these  particular  solutions.' 

[1981]  has  attempted  to  extend  D’Yakonov’s  solu¬ 
tion  to  evaluate  the  fields  at  the  surface  of  the  earth 
using  the  circular  cylinder  as  a  target.  However,  in 
his  assumed  form  of  the  scattered  fields  he  has  intro¬ 
duced  a  symmetry  in  the  angle  0  about  <t>o  that  is  not 
actually  present  (see  his  equation  (21)  and  his  Figure 
2).  Even  if  this  were  corrected  by  introducing  the 
terms  sin  n(0  -  <t>o),  the  result  would  still  not  be  cor¬ 
rect  except  for  the  observation  point  directly  above 
the  buried  cylinder.  Also,  as  we  have  already  ob¬ 
served,  it  is  not  always  essential  when  considering  the 
scattering  problem  to  include  the  half  space  problem. 
Indeed,  much  valuable  insight  can  be  achieved  in  an 
economic  manner  if  it  is  not  included.  Modal  solu¬ 
tions  for  scattering  from  lossy  circular  cylinders  in 
lossy  half  space  excited  by  a  line  source  are  relatively 
straightforward  to  evaluate  numerically  using 
modern  digital  machines.  However,  the  numerical 
computations  become  more  complex  when  a  finite 
antenna  system  (neither  infinite  nor  infinitesimal)  is 
introduced  at  the  surface  [Burre/f  and  Pthrs,  1979]. 
There  have  been  many  analyses  involving  infinitesi¬ 
mal  dipole  sources  located  at  the  surface  of  a  half 
space.  Many  of  these  include  reflections  from  multi¬ 
ple  layers  in  the  lossy  half  space  [IVai/,  1958].  There 
have  been  few  numerical  solutions  involving  finitc- 
sized  antennas,  particularly  for  the  current  distri¬ 
bution  and  the  impedance,  such  as  a  resonant  dipole 
at  the  surface  of  the  lossy  half  space,  primarily  be¬ 
cause  of  computation  costs.  This  has  been  approxi¬ 
mated  by  Davis  [1979].  There  is.  however,  a  conti¬ 
nued  effort  to  generate  new  techniques  which  will 


ultimately  eliminate  this  difficulty.  Howard  [1975], 
for  example,  uses  fast  Fourier  transform  (FFT)  rou¬ 
tines  to  obtain  the  scattered  fields  of  a  small  radius 
conducting  cylinder  immersed  in  a  lossy  half  space. 
There  is  at  least  one  other  modal  solution  that 
should  be  included  in  this  discussion. 

Wait’s  [1959]  treatment  of  the  scattered  fields 
from  an  imperfectly  conducing  dike  should  also  be 
included  as  a  modal-based  solution.  His  approach  is 
a  rather  elegant  preface  to  concepts  now  being  intro¬ 
duced  in  the  geometrical  theory  of  diffraction  (GTD) 
for  the  treatment  of  thin  dielectric  slabs  and  absorb¬ 
ers  placed  on  conducting  surfaces.  We  will  later 
update  this  technique,  show  some  newer,  remarkably 
accurate  results,  and  propose  some  possible  exten¬ 
sions  to  this  approach.  The  modal  solution  is  appli¬ 
cable  primarily  to  cylindrical  shapes.  With  current 
computer  technology  there  is  no  serious  limitation 
on  the  size  of  the  cylinder  that  can  be  treated.  In  this 
solution  the  incident  wave  is  expanded  in  a  set  of 
cylindrical  modes,  and  boundary  conditions  are  ap¬ 
plied  to  each  mode  individually  to  evaluate  the  coef¬ 
ficients  of  the  scattered  cylindrical  modes.  Summing 
the  fields  of  each  mode  then  gives  the  desired  scat¬ 
tered  fields. 

Our  recent  computations  have  involved  a  number 
of  solutions  for  the  SAF  of  circular  cylinders.  Some 
examples  arc  shown  in  Figures  2-4.  Figure  2  gives 
the  SAF  for  1-m-radius  tunnels  of  different  depths 
and  ambient  media.  It  is  left  as  an  exercise  to  check 
the  guidelines  for  depth  and  the  media  electrical  pa¬ 
rameter  dependence  given  earlier.  Figure  3  gives 
SAF’s  for  lossy  cylinders  immersed  in  a  lossy 
medium.  The.se  might  represent  tunnels  filled  with 
water  or  debris.  Figures  4  and  5  give  the  SAP's  for 


.1.  Scatlering  altenuaiion-  funoiion  for  lossy  dielectric  cyl¬ 
inders  in  a  lossy  medium  for  parallel  electric  polariralion. 
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Fig.  4.  Scallering  attenuation  function  for  air-Olled  cylinder 
with  a  perfectly  conducting  wire  at  the  tunnel  floor. 

air-filled  cylinders  with  a  wire.  The  wires  might  rep¬ 
resent  a  rail  or  a  power  line  in  the  tunnel.  Let  £| 
denote  the  field  of  the  line  source  near  the  dielectric 
cylinder  without  wire.  Now  let  fij  denote  the  field  of 
a  1-A  line  source  located  within  the  dielectric  cylin¬ 
der  at  the  position  where  we  will  later  place  the  wire. 
The  solution  for  makes  use  of  the  addition  theor¬ 
em  [Harrington,  1961].  By  superposition,  the  field  of 
the  exterior  line  source,  radiating  in  the  presence  of 
the  dielectric  cylinder  with  interior  wire,  is  given  by 
£  =  £i  +  lEj.  The  current  I  induced  on  the  thin 
wire  is  determined  by  forcing  £  to  vanish  on  the 
surface  of  the  wire.  The  current  is  assumed  to  be 
distributed  uniformly  around  the  circumference  of 
the  wire.  Finally,  Figure  6  shows  the  SAF  for  a  con- 


Fig.  6.  Scattering  attenuation  function  for  perfectly  conducting 
cylinder  and  perpendicular  polarization  in  relatively  loss  free 
ground. 

ducting  cylinder  immersed  in  a  nearly  lossless 
medium.  The  ripple  in  this  case  is  caused  by  the 
creeping  wave  which  is  almost  negligible  even  for  this 
very  lossless  medium. 

ALTERNATE  TECHNIQUES  AVAILABLE  FOR 
SCATTERING  FROM  PENETRABLE 
CYLINDERS 

There  are  perhaps  three  additional  techniques 
available  to  evaluate  scattered  fields  from  penetrable 
cylinders.  These  include  (1)  moment  method  or  inte¬ 
gral  equation  solutions,  (2)  optical  solutions,  and  (3) 
hybrid  solutions.  Each  solution  has  its  applications 
and  its  limitations,  but  they  have  not  been  developed 
at  this  time  to  their  ultimate  capacity  for  targets  im¬ 
mersed  in  a  lossy  medium. 


Im  RADIUS  CYLINDER: 


Fig.  3  Scallering  ailenualion  function  for  air-filled  cylinder 
with  a  perfectly  conducting  wire  at  the  roof  of  the  tunnel. 


INTEGRAL  EQUATION  SOLUTIONS 

There  are  two  basic  approaches  using  integral 
equations  for  the  scattered  fields  of  penetrable  ob¬ 
jects.  One  of  these  uses  surface  equivalence  currents 
(J  =  nxH  and  M  =  E  x  »?)  [Harrington,  1961]  and 
requires  the  solution  of  coupled  integral  equations, 
one  for  the  ambient  medium  and  the  other  for  the 
penetrable  object  [Andreason,  1965].  Matching  the 
fields  at  the  surface  allows  the  unknown  coefficients 
to  be  evaluated.  This  is  quite  similar  to  the  Sommer- 
fcld  analysis  of  the  h.ilf  space.  The  other  approach 
uses  the  so-called  vi  vtc  polarization  currents  to  be 
discussed  later  [Rhodes,  1953;  Richmond,  1965, 
1966].  One  of  the  major  distinctions  is  that  in  the 
first  case  the  boundary  conditions  are  satisfied  at  the 
surface  and  in  the  second  case  they  arc  satisfied 
throughout  the  volume.  The  choice  of  solution  would 
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be  associated  largely  with  the  electrical  parameters, 
size,  and  shape  of  the  penetrable  body,  and  the  de- 

v/«zulU  VfT  *in  tjic  matrix  t#>  he 

inverted. 

For  example,  if  a  thin  layer  is  to  be  treated,  then 
the  polarization  current  approach  would  be  most  ap¬ 
propriate.  Alternately,  if  a  large  perturbation  at  an 
interface  is  to  be  analyzed,  then  the  surface  current 
would  be  more  appropriate.  Smaller  objects  or  ob¬ 
jects  with  a  small  dimension  in  terms  of  wavelength 
would  require  a  smaller  matrix  if  the  volume  polari¬ 
zation  current  approach  is  used.  Since  most  moment 
method  solutions  are  considered  low-frequency  tech¬ 
niques,  then  it  seems  appropriate  to  place  emphasis 
on  the  polarization  current  density  approach.  An  ex¬ 
ception  to  this  might  be  a  gentle  perturbation  or 
‘swell’  on  an  interface. 

There  is  also  an  integral  equation  solution  that  has 
been  designated  as  the  T  matrix  used  originally  by 
Waterman  [1965]  for  the  analysis  of  electromagnetic 
scatterers.  These  solutions  incorporated  an  extended 
boundary  condition  so  that  external  to  this  bound¬ 
ary,  outgoing  cylindrical  (or  spherical)  waves  could 
be  used.  Mei  et  al.  [1979]  used  a  similar  solution 
except  that  they  introduced  a  (ini!  -  Jitfercncc  solu¬ 
tion  to  hnd  the  trial  functions  internal  to  the  ex¬ 
tended  boundary.  Again,  these  types  of  solutions 
would  minimize  the  number  of  match  points  where 
the  scatterer  better  fits  the  circular  or  spherical  ex¬ 
tended  boundary. 

Next  we  will  outline  the  solution  for  a  lossy  dielec¬ 
tric  of  rather  general  cross  section  using  the  polariza¬ 
tion  current  density  approach. 

MOMENT  METHOD  l  OR  ELECTRIC  LINE  SOURCE 
NEAR  A  LOSSY  DIELECTRIC  CYLINDER 
USING  THE  POLARIZATION  CURRENT 
APPROACH 

Consider  an  inlinite  electric  line  source  with  a  uni¬ 
form  time  harmonic  current  /q  ■  When  radiating  in  a 
homogeneous  medium  with  complex  parameters  ft 
and  e,  this  source  generates  the  following  incident 
field: 

E' *  -ifo»n'A'„()7),)/2a  (4) 

where  ii  and  v  denote  the  intrinsic  impedance  and 
propagation  constants  of  the  ambient  medium  and 
Ko  denotes  the  modified  Bessel  function.  I'he  line 
source  is  UKuted  at  (.v,,  y,)  and  is  parallel  with  the  : 
axis.  The  observer  is  located  at  (.v,  >j.  and 

p,  -  [(X-.\,)'  +()•->■,)']'  ' 


Now  consider  the  electric  line  source  radiating  in  the 
presence  of  a  nearby  parallel  dielectric  cylinder.  The 
field  in  this  case  is  denoted  by  E.  where 

E  =  £'  +  (6) 

Since  the  field  is  ;-pularized,  scalar  notation  is  suit¬ 
able.  The  lield  £•/'  scattered  by  the  dielectric  cylinder 
may  be  generated  by  an  equivalent  polarization  cur¬ 
rent  density  given  by 


J ^  Mf  ,  -  i  )i--  (7) 

radiating  in  the  homogeneous  ambient  medium  (/ic). 
The  dielectric  cylinder  has  parameters  /i,  and  t,,  and 
it  may  be  inhomogeneous,  bui  we  assume  /i,  =  /j. 
The  current  7,,  exists  only  in  the  dielectric  region 
with  permittivity  Since  the  field  E  is  unknown,  the 
current,  is  also  unknown. 

The  scattered  field  is  given  by 


E%x.  .»■)  = 


2n 


y'}K„(yp)  dx‘  dy’ 


(8) 


where  the  integration  extends  over  the  cross- 
sectional  area  of  the  cylinder  and 


p  =  [(x-\')-  4  (9) 


From  (6)  (8)  we  can  derive  the  following  integral 
equation ; 


J{x'.  ,v')Ko(;7>)  dx'  dy’  =  E\x.  >j  (10) 

where  d  =  J,, .  Now  let  us  expand  d  in  a  series  of 
basis  functions  as  follows: 


d(.x..v)  ^  nr  rr 
doXci-f.)  2n  JJ 


d(x.  ,v) --- V  >•)  (II) 

From  (10)  and  (I  I), 


F.<x.  y) 


SIJ 


E,(  x',  >’■)£„(•, 7>)  dx'  dy 


>'] 


=  FAX.  y) 


(12) 


To  determine  the  constants  let  us  multiply  (12) 
by  a  weighting  function  (;„(\.  \  |  and  integrate  over 
the  cross-sectional  area  of  the  cylinder.  In  this 
manner  we  obtain  the  following  system  of  simulta¬ 
neous  linear  equations: 


m-1.  2.  N 

N  I 


<*1.1  X.  >  lE,(  v.  l  )  «/v 
jiiitr.,  -  c) 


ny 

2n 


.»•) 


(13) 


(5) 


•  fjlx'.  >•■)£„(', 7>)  ds  ds 


(14) 


y.  =  jj  GJx,  y)E\x,  y)  ds 


Numerical  integration  techniques  are  generally  re¬ 
quired  to  evaluate  the  integrals  in  (14)  and  (IS). 
Matrix  inversion  yields  the  solution  to  (13),  and  then 
the  scattered  held  is  calculated  via  (8)  and  (11). 

MOMENT  METHOD  FOR  MAGNETIC  LINE  SOURCE 
NEAR  A  LOSSY  DIELECTRIC  CYLINDER 

Consider  an  infinite  magnetic  line  source  with  uni¬ 
form  time  harmonic  current  distribution  Iq.  When 
radiating  in  a  homogeneous  medium  with  parameters 
H  and  £,  this  source  generates  the  following  incident 
field; 

2nri 

tidH,  y/o(y-y.)K,(y#),) 

K - ^  -  (17) 

y  2itp, 

nSH,  y/o(x-x.)K,(yp.) 

y  cx  2np, 

The  field  scattered  by  the  parallel  dielectric  cylinder 
(or  radiated  by  the  polarization  current  J  in  a  homo¬ 
geneous  ambient  medium)  is  given  by 

HK*.  >')  =  I;  II  CO’-/)-'.!*',  y') 

KAyp) 

-{x-xV^x-,y')}-^dx-  dy'  (19) 

P 


J  =  ytu(£,  -r)F. 

From  (20)  and  Maxwell's  equations. 


J. 


(t.-£)  SH, 
£,  dy 


—  j(u(Ci  c)Ey  ■ 


From  (l9)-(22). 


i(yF)  . ,  , , 

- dx  dy 

P 


The  incident,  scattered,  and  total  fields  are  related  as 
follows: 

H  -  H'  +  H‘  (24) 
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From  (23)  and  (24)  we  obtain  the  following  integral 
equation; 

ax  J  p 

Now  let  us  expand  the  unknown  function  //  in  a 
series  as  follows; 

H(x,  y)  =  H.(x,  >•)  =  X  C.  F,{x,  y)  (26) 

where  the  functions  F,  are  called  the  basis  functions. 
From  (25)  and  (26), 

z  c.  f.(x,  y)  =  «‘(x,  >•)  +  Xc.  II 

r  dF,  dF,l  /C,(yp) 

•  I  O'  -  y')  ^  +  (X  -  x  )  dx'  dy  (27) 

To  determine  the  constants  C„.  we  multiply  (27)  by  a 
weighting  function  CJx,  y)  and  integrate  over  the 
cross-sectional  area  of  the  dielectric  cylinder.  This 
yields  the  following  system  of  simultaneous  linear 
equations: 

N 

lc.z„=v„  m=1.2.-N  (28) 

K*  I 

z-  -  IJcjx.  JX.  y)  *-  i  J|  IJ  ^  <rjx.  ,1 

•  [(y  -  y')  ^  +  (X  -  X')  ^1  ds'  ds  (29) 

L  ^y  ^x  J  p 

»''-  =  ||cjx,>)W'(.v,>)ds  (30) 

Equation  (28)  is  solved  by  matrix  inversion  or  equiv¬ 
alent  techniques,  and  then  the  scattered  field  is  calcu¬ 
lated  via  (23)  and  (26). 

Figure  7  shows  the  scattering  attenuation  function 
as  defined  in  (I)  for  a  square  cylinder  as  obtained 
with  the  general  approach  just  outlined.  Hohman 
also  made  some  scattered  field  computations  from  a 
square  conducting  medium,  but  these  did  not  extend 
into  the  resonance  region. 

Use  of  the  SAF  of  Figure  1  in  the  model  of  Figure 
1  has  yielded  reasonable  agreement  with  measured 
results  as  obtained  using  a  video  pulse  radar  [Davis, 
1979], 

INSULATED  FINITE  LENGTH  WIRES 
IN  A  CONDUCTING  MEDIUM 

Conducting  finite  length  wires  immersed  in  a  con¬ 
ducting  medium  can  be  readily  analyzed  using  con- 
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Fig.  7.  Composite  scattering  function  of  2-m-square  air-filled 
tunnel  obtained  by  a  polarization  current  moment  method  solu¬ 
tion.  fc'i  -  £||  represents  (he  SAF  for  a  crossed  dipole  antenna  (see 
equation  (3)). 

ventional  moment  method  techniques.  One  such 
result  has  already  been  noted.  We  have  also  con¬ 
sidered  infinite  wires  immersed  in  a  circular  tunnel 
using  modal  techniques.  The  solution  to  be  discussed 
here  is  that  of  an  insulated  finite  length  wire.  A 
coaxial-like  mode  is  assumed  where  the  volume  po¬ 
larization  currents  are  placed  in  the  dielectric  shell 
and  surface  currents  are  on  the  surface  of  the  wire. 
The  electric  field  in  the  dielectric  shell  is  assumed  to 
be  radial  yielding  polarization  currents 

•/,  =  -  (£2  -  c)/'(0  (31) 

2nr.i  p 

where  /(/)  is  the  current  on  the  wire.  Designating  this 
current  in  the  dielectric  shell  as  a  tubular  expansion 
dipole,  one  needs  to  include  the  mutual  impedance 
between  the  nth  tubular  exp.snsion  dipole  and  the 
mth  filamentary  dipole  (in  the  lossy  medium)  in  the 
mutual  impedance  Z„,.  The  moment  method  solu¬ 
tion  then  proceeds  in  a  normal  fashion.  The  reader  is 
referred  to  Richmond  and  Newman  [1976]  for  further 
details.  This  approach  has  been  tested  by  evaluating 
the  impedance  of  insulated  wire  antennas.  A  typical 
example  is  shown  in  Figures  8  and  9.  This  same  ap¬ 
proach  has  been  used  to  find  the  scattered  fields  from 
an  insulated  conducting  wire  in  a  homogeneous 
medium.  In  this  case,  we  present  the  scattering  at¬ 
tenuation  function  results,  as  incorporated  in  the  re¬ 
ceived  voltage  at  the  terminals  of  the  receiving  an¬ 
tenna  shown  in  Figure  1.  The  source  in  this  case  is  a 
Gaussian  pulse  which  contains  a  rather  wide  spec¬ 


trum.  The  steps  in  the  analysis  include  taking  the 
FFT  of  the  transmitting  pulse,  evaluating  the  compo¬ 
nent  terms  illustrated  in  Figure  I,  and  taking  the 
inverse  FFT  to  obtain  the  received  voltage  wave¬ 
forms  now  shown  in  Figure  10.  The  spiectrum  of  sev¬ 
eral  of  these  waveforms  is  shown  in  Figure  1 1,  where 
‘tunnel  length'  represents  the  length  of  the  insulated 
buried  wire. 

HYBRID  SOLUTIONS 

Hybrid  solutions  have  been  used  to  analyze 
various  geometries  that  may  not  be  treated  by  the 
usual  moment  method  solutions.  They  have  been 
used  to  treat  the  infinite  conducting  wedge  [Burnside 
et  ai,  197S]  and  an  initial  solution  for  an  infinite 
dielectric  wedge.  They  have  not  yet  been  used  for  the 
treatment  of  objects  immersed  in  a  lossy  medium. 

The  hybrid  solution  is  a  combination  of  a  moment 
method  solution  with  another  solution  where  the  un¬ 
known  fields  out.sidc  a  restricted  region  can  be  rep¬ 
resented  by  a  field  of  known  functional  form  multi¬ 
plied  by  an  unknown  constant.  One  of  the  original 
solutions  took  the  form  illustrated  in  Figure  12.  The 
currents  on  the  surface  of  the  conducting  wedge  out¬ 
side  the  encircled  region  can  be  written  in  the  form 

Ce 

J“Ti>  =  J'  +  J'  +  (a  X  ■)  — —  (32) 

P 

where  J'  -  n  x  If  and  =  >5  x  H'.  H'  and  H'  are 
the  incident  and  reflected  magnetic  field  intensities 
respectively;  the  only  unknown  is  the  constant  C. 


Fig.  8.  The  conductance  of  an  insulated  dipole  plotted  versus  its 
length  in  wavelengths. 
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The  surface  currents  inside  the  encircled  region  are 
unknown  and  are  evaluated  along  with  the  constant 
C  using  conventional  moment  method  techniques. 
Harrington  [1968]  gives  as  an  integral  equation  for 
the  total  current  density  J  : 

J  -  -  +  z  .  V  X  ^  j  JH!j‘'(k|p  -  p'\)dl'^  (33) 

where  L  is  taken  about  the  perimeter  of  the  conduc¬ 
ting  surface.  Inside  the  encircled  region  the  currents 
can  be  represented  using  subsectional  bases.  Outside 
the  encircled  region  they  are  represented  by  the  GTD 
currents  given  by  (32).  For  the  lossless  media,  this 
will  involve  integrations  along  the  surface  of  the  form 

ik 

-  —  H\^\k\p-p'\)dp'  (34) 

fit  ^ 

where  p'  is  taken  along  the  surface  of  the  wedge. 

As  Burnside  et  al.  [1975]  observe,  this  integration 
need  not  be  carried  to  oo  because  of  the  decaying 
nature  of  the  Hankel  function. 

Next,  one  can  consider  the  case  of  the  conducting 
wedge  in  a  lossy  medium  as  a  first  step  in  considering 
lossy  wedges  immersed  in  a  lossy  medium  as  the  ex¬ 
ternal  medium  becomes  lossy  (/k-»  —a—  jk)\  this  in¬ 
tegral  becomes  negligible  for  much  smaller  values  of 
p,  and  thus  the  solution  is  simplified.  If  the  loss  is 
more  than  marginal,  then  indeed  one  could  find  the 
fields  diffracted  by  the  edge  of  the  wedge  without 
using  the  hybrid  approach,  unless  one  is  particularly 
interested  in  the  fields  along  the  surface  of  the  wedge. 
Aas  [1975]  has  used  additional  terms  in  of  the 
form  and  Cje'^^lpi  to  include  the  cases  of 


DIPOLE  LENGTH  L/X 


Fig.  9.  The  tiMCcpIance  of  (he  iniulmed  dipole,  shown  in  Figure 
8,  piolled  versus  its  length  in  wavelengths. 


Fig.  to.  Signals  received  on  a  50-m  orthogonal  dipole  antenna 
system  SO  m  from  300  m  long,  O.l-m-radius  tunnel  containing  a 
0.01-m-radius  wire  on  its  axis,  (a  and  b)  The  early  lime  behavior 
and  the  complete  received  signal  respectively  for  a  O.S-fis  Gaussian 
input  pulse,  {c-d  and  e-f)  The  signals  for  t-;<s  and  2-;iS  input 
pulses  respectively. 

grazing  incidence  and  also  polarization  parallel  to 
the  surface. 

Wu  and  Tsai  [1977]  have  adapted  this  approach 
for  the  case  of  a  dielectric  wedge.  They  have  used 

JOTO  =  J<  +  jr  (35) 

They  have  included  only  first-order  geometrical 
optics  fields  for  the  region  far  away  from  the  tip  of 
the  wedge,  primarily  because  there  is  not  yet  a  solu¬ 
tion  for  the  diffracted  fields  of  a  dielectric  wedge. 

At  first  glance,  it  would  appear  that  the  coupled 
integral  equation  approach  would  be  more  suitable 
for  the  hybrid  approach  since  it  would  be  more  suit¬ 
able  for  defining  a  wave  on  the  surface  of  the  wedge 
comparable  to  the  hybrid  analysis  of  the  conducting 
wedge.  They  observe  that  the  solution  could  be  im¬ 
proved  by  including  multiple  reflections  internal  to 
the  wedge.  They  also  note  that  it  would  be  improved 
when  the  tip  diffraction  term,  once  its  form  is  known, 
is  included. 

OPTICAL  .SOLUTIONS 

The  optical  solutions  generally  analyze  scattering 
using  ray  techniques  such  as  conservation  of  energy 
in  a  ray  tube,  plane  wave  refraction  and  reflection 
coefficients,  caustics,  etc. 

MODIKIHD  GEOMETRICAL  OPTICS  METHOD 

For  example,  some  of  the  earliest  solutions  of  this 
type  focused  attention  on  the  scattered  fields  of  di- 
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Fig.  II.  Magnitude  of  transmission  transfer  function  for  50-m- 
long.  0.002-m-radius  uninsulated  dipole  antenna  system  SO  m 
from  a  finite  length  tunnel  for  various  tunnel  lengths.  The  tunnel 
has  a  radius  of  0.1  m  and  contains  a  0.01-m-radius  wire  at  its 
center. 


electric  spheres  and  cylinders  [Kouyoumjian  et  al, 
1963;  Peters  et  at.,  1965].  Analyses  of  this  type  were 
designated  as  the  modihed  geometrical  optics 
method.  Some  of  the  ray  mechanisms  that  were  con¬ 
sidered  are  illustrated  in  Figure  13.  Ray  1  in  Figure 
13  is  an  externally  reflected  ray,  ray  2  is  the  axial 
internally  reflected  ray,  and  ray  3  represents  two 
nonaxial  rays.  One  of  these  is  shown  in  Figure  14 
and  is  known  as  the  glory  ray.  It  has  an  apparent 
focus  or  caustic  located  a  distance  pi  from  the  sur¬ 
face  as  shown.  Figure  15  shows  the  geometry  of  the 
stationary  or  rainbow  ray.  The  analysis  included  the 


F  ig.  12.  Geometry  for  hybrid  lolulion  for  (he  conducting  wedge. 


electrical  path  lengths,  the  caustics,  and  multiple  in¬ 
ternally  reflected  rays.  Figure  16  shows  the  compari¬ 
son  of  the  approximate  echo  width  of  several  cylin¬ 
ders  whose  relative  permittivities  are  0.50  and  0.75 
with  those  obtained  from  the  exact  solution.  These 
examples  were  chosen  since  they  might  be  more  rep¬ 
resentative  of  a  sewer  pipe  or  a  plastic  gas  pipe, 
tunnel,  etc.,  than  other  values  that  have  also  been 
obtained.  The  analysis  becomes  more  complex  as 
e,  =  Ci/cq  shown  in  Figure  13  becomes  large  with 
respect  to  1.  landa  and  Planus  [1969]  have  criticized 
this  approach  and  claimed  that  surface  waves  are 
needed  to  complete  the  picture.  There  is  one  other 
mechanism  that  is  possible,  a  lateral  (or  up.  over,  and 
down)  wave.  This  would  be  applicable  when  the 
index  of  refraction  for  the  cylinder  is  greater  than 
that  of  the  ambient  medium.  At  any  rate  for  e,  = 
2.56,  Barrick  [1968]  included  the  following  rays:  (1) 
externally  reflected  ray,  (2)  single-bounce  axial  ray, 
(3)  triple-bounce  axial  ray,  (4)  single-bounce  glory 
ray,  and  (5)  single-bounce  stationary  ray.  He  ob¬ 
tained  reasonably  good  agreement  between  the  mod¬ 
ified  geometrical  optics  solution  and  the  exact  solu¬ 
tion  except  for  raggedness  in  the  frequency  depen¬ 
dent  curves  of  the  exact  solution.  He  reported  that 


I'ig.  14.  Wive  fruni  of  single-bounce  glory  riy. 
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Fig.  ISa.  Equiphase  surfaces  associated  with  single-bounce 
stationary  ray. 


Fig.  1 5b.  Enlarged  figure  about  point  \  for  the  equiphase  sur¬ 
face  of  single-bounce  stationary  ray. 


the  addition  of  small  losses  eliminated  this  ragged 
characteristic.  Including  other  rays  of  the  type  listed 
did  not  introduce  such  a  behavior,  and  he  attributed 
it  to  a  much  higher  order  periodic  type  of  trapped 
mode.  This  is  probably  the  lateral  wave  mechanism. 

This  type  of  solution  has  been  applied  to  a  variety 
of  dielectric  materials  where  the  targets  are  generally 
circular  cylinders  and  spheres.  There  is,  however, 
nothing  that  restricts  it  to  these  geometries.  It  could 
be  applied  for  almost  any  shape  exclusive  of  caustic 
regions  where  special  techniques  are  needed  to  evalu¬ 
ate  the  fields.  While  all  analyses  thus  far  have  con¬ 
sidered  only  lossless  targets  in  a  lossless  medium 
losses  can  be  accounted  for  in  terms  of  an  attenu¬ 
ation  along  the  pertinent  ray  paths,  etc.  These  ray 
paths  may  lie  in  the  ambient  medium  and  on  or 
inside  the  scatterer.  The  introduction  of  loss  may  also 
generate  additional  ray  paths.  Of  course,  rellection 
coefficients,  etc.,  will  be  modified.  It  should  be  noted 
that  such  solutions  would  generally  be  simpler  since 
multiple  internal  reflections  in  a  lossy  cylinder  would 
be  attenuated  and  would  reduce  the  number  of  inter¬ 
nal  rays  that  need  be  included  in  the  solution. 

GEOMETRICAL  THEORY  OF  DIFFRACTION 
SOLUTION 

A  more  recent  optical  solution  makes  use  of  some 
of  the  concepts  of  edge  dilTraction  [Burnside  and 
Pathak,  1981].  Consider  the  diffracted  fields  from  the 
edge  of  a  conducting  half  plane  given  by 

W\  jU>\ 

or  I  «  I  or  l[D(.i,  .1',  il>  -(p‘)±  D(.v,  s',  <t>  +  0'H  (36) 

where  £'  and  H'  are  the  fields  of  the  incident  wave.  D 
is  the  diffraction  coeflicient,  and  the  geometry  is 


given  in  Figure  17.  The  incident  fields  are  shielded  at 
<l>  —  <f>'  =  n,  which  is  called  the  shadow  boundary. 
The  image  (or  reflected  fields  from  the  half  plane) 
fields  do  not  exist  for  (j>  +  (j>'  =  n.  Thus  a  reflection 
boundary  exists  at  <f)  -(-</)'  =  n.  The  diffraction  coef¬ 
ficients  D(s,  s',  -  Ip')  “od  D[s,  s.  Ip  +  Ip')  are  inti¬ 

mately  associated  with  the  incident  and  reflection 
.shadow  boundaries  respectively.  The  diffracted  fields 
given  by  f)(s.  s',  ip  +  ip')  are  proportional  to  the  dis¬ 
continuity  of  the  fields  at  the  shadow  boundary  (SB), 
and  the  diffracted  fields  given  by  D{s,  s',  ip  -f  iP')  arc 
proportional  to  the  discontinuity  at  the  reflection 
boundary  (RB). 


•  •  •  •  AP''ROXIUATE 

- EXACT  I  PARALLEL  AND  PERPENDICULAR 


Fig.  Ifi.  r»mpAri<H<n  of  exact  anti  apprtiximale  backxcalieT 
echo  width  of  an  infinite  circular  dielectric  cylimler  tor  rctalivc 
permittivity  r.,  «  0.50, 0.75. 
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Fig.  17.  Edge-dilTracted  ray  geometry.  RB  is  reflection  boundary; 

SB  is  shadow  boundary. 

Suppose  now  we  permit  (he  half  plane  to  be  thin 
and  penetrable.  Then  there  is  a  partial  reflection  and 
a  partial  transmission  at  the  half  plane  as  represented 
by  Figure  18.  Now  (36)  takes  the  form 

I  [(1  -  T)iJj(.s,  s',  <t>  -  4>') 

+  R.D(s.  s',  0  +  <t>'n  (37) 


halt  planc 

Fig.  IS.  Optical  houndarin  for  dielectric  half  plane 


k 


Fig.  19.  Geometry  for  CslO  moment  melhoJ  (MM)  com* 
pa  risen. 


where  (!  —  T)  =  R  represents  the  discontinuity  of 
fields  at  the  shadow  and  reflection  boundaries  respec¬ 
tively.  This  type  of  solution  has  been  applied  to  the 
case  of  a  thin  dielectric  slab.  Burnside  and  Pathak 
[1981]  have  calculated  the  far  field  radiated  fields  for 
a  number  of  two-dimensional  geometries.  An  exam¬ 
ple  of  this  geometry  is  given  in  Figure  19.  A  typical 
result  is  shown  in  Figure  20  where  this  GTD  solution 
is  compared  with  the  moment  method  solution. 

We  observe  that  the  results  obtained  herein  do 
confirm  the  accuracy  of  the  method.  The  fact  that 
there  is  agreement  over  the  entire  pattern  also  re¬ 
moves  the  potential  restriction  discussed  by  Wait 
[1959]  concerning  the  location  of  the  source  and  ob¬ 
servation  point.  This  is  accomplished  by  the  appro¬ 
priate  modification  of  the  diffracted  fields  associated 
with  both  the  shadow  and  the  reflection  boundaries. 

A  significant  factor  here  is  that  the  scattered  fields 
can  indeed  be  represented  by  a  discrete  shadow 
boundary  and  by  the  same  form  of  diffracted  fields  as 
for  the  conducting  shape.  If  the  same  basic  postulate 
is  made  for  penetrable  wedge  geometries  with  differ- 


Fig.  20.  GTD-MM  o>n)p;iri«>ii  wilh  clix'lric  line  wurce  using 
geomtiry  shown  in  Figure  !•>  wilh  I.  -  20/.  f)  *  0.05/,  p,  2  02 
(Lefll  0,-20  ami  irighll  0,  40. 
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wedge. 


ent  angles,  it  may  become  possible  to  construct  solu¬ 
tions  for  more  general  shapes.  In  the  following  para¬ 
graphs,  we  outline  some  of  the  considerations  that 
might  be  used  to  obtain  GTD  type  solutions  to  more 
complex  geometries.  The  reader  is  cautioned  that 
these  are  merely  presented  here  us  untested  postu¬ 
lates. 

The  principle  limitation  on  the  results  shown  in 
Figure  20  is  the  thickness  of  the  slab.  For  thicker 
slabs  or  for  subsurface  geometries  that  might  be  of 
interest,  the  index  of  refraction  and  the  angle  of  inci¬ 
dence  will  introduce  additional  limitations.  Consider, 
for  example,  the  geometry  of  Figure  21.  For  this  case, 
it  might  appear  that  the  boundary  conditions  are 


'^(SBIFOR 

FAST 

WAVE 


Fig.  23.  Phase  velocilies  of  a  wave  along  the  surface  of  the  wedge 
of  Figure  22. 


essentially  those  that  have  already  been  stated  and 
used  in  (37),  and  indeed  for  diffracted  fields  in  the 
vicinity  of  the  reflection  boundary  (where  the  second 
term  of  (37)  is  dominant)  this  may  represent  a  valid 
approximation.  However,  for  fields  in  the  vicinity  of 
the  shadow  boundary  (where  the  first  term  in  (37)  is 
dominant)  this  approach  is  clearly  erroneous.  Also, 
since  the  phase  velocities  of  the  ‘ray  optical  fields'  in 
the  two  media  must  differ,  it  is  impossible  to  satisfy 
the  boundary  conditions  along  the  surface  by  such 
fields.  Thus  we  now  must  introduce  the  additional 
complication  of  some  type  of  'trapped'  or  lateral 
wave  field.  Thus  the  possibility  of  a  Zenncck  type 
wave  [Hi//  and  IVait,  1978,  1980]  must  be  con¬ 
sidered.  These  fields  would  be  represented  by  ‘com¬ 
plex  rays.'  This  subject  has  been  studied  in  consider¬ 
able  depth  by  Felson  and  his  co-workers  [Fc/son, 
1981]. 


Fig.  24.  A  penelrablc  wedge  geometry  where  the  GTD  solution 
should  give  good  results. 


Fig.  22.  Grexinf  incideitce  on  the  face  of «  pcnetriMe  wedge. 
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Fig.  2S.  Suggested  ray  configuration  for  a  simplistic  fault  geom¬ 
etry  where  current  GTD  solution  is  appropriate. 

Looking  further  at  the  geometry  of  the  vertical 
interface  of  Figure  21  as  illustrated  in  Figure  22,  we 
observe  that  the  phase  velocity  of  the  wave  at  the 
surface  on  the  low  index  of  refraction  side  has  the 
appearance  of  a  fast  wave.  Its  phase  velocity  must 
have  a  component,  not  only  along  the  vertical  inter* 
face  but  also  perpendicular  to  it.  It  would  conse¬ 
quently  appear  as  a  shadow  boundary  as  shown  in 
Figure  23. 

The  launching  properties  of  such  a  field  should  be 
studied  carefully  using  the  hybrid  solutions  already 
discussed.  Once  these  are  adequately  understood, 
they  would  be  used  as  a  tool  for  studying  more  com¬ 
plex  geometries. 

Note  also  that  the  various  boundaries  are  associ¬ 
ated  with  different  wedge  surfaces.  For  example, 
(SB)]  is  clearly  associated  with  the  vertical  sufacc  of 
Figure  21.  Thus  once  the  equivalence  currents  are 
evaluated  for  the  wedge,  one  could  assume  that  the 
equivalence  currents  on  the  vertical  surface  would 
contribute  to  the  fields  (SB)]  and  its  associated  dif¬ 
fracted  fields.  One  could  further  refine  this  model  by 
neglecting  the  optical  fields  on  the  vertical  edge  and 
repeating  the  hybrid  solution.  In  this  manner,  one 
could  hope  to  construct  a  solution  for  a  dielectric 
wedge  in  a  GTD  format. 

In  Figure  21  we  recall  that  there  now  are  two 
shadow  boundaries  at  different  angles.  From  a  geom¬ 


etric  optics  viewpoint,  there  is  a  void  between  these 
two  boundaries.  This  is  clearly  a  more  complex  case. 
It  becomes  even  more  complex  if  the  line  source  is 
not  in  the  far  field  since  the  fields  of  rays  associated 
with  (SB)}  are  not  homogeneous  because  they  are 
dependent  on  various  incidence  angles  at  the  surface. 
One  could  make  various  postulates  as  to  the  form  of 
these  diffracted  fields.  However,  at  this  time,  not 
enough  has  been  achieved  to  be  realistic.  Perhaps  the 
proper  approach,  barring  the  development  of  an  ade¬ 
quate  modal  solution,  would  be  to  extend  the  hybrid 
solution  and  seek  GTD-like  repre.sentations  with  the 
appropriate  uniform  behavior  to  represent  them  as 
rays. 

Figure  24  illustrates  an  example  where  reasonable 
results  might  be  obtained  for  the  space  designated  as 
If,  using  the  approach  given  by  Burnside.  The  overly 
stringent  requirement  that  i;,  <  if}  implies  that  there 
are  no  multiply  reflected  rays.  The  diffracted  fields 
from  any  internal  shadow  boundaries  would  tend  to 
remain  inside  the  wedge.  The  external  diffracted 
fields  would  be  those  of  the  conducting  wedge  multi¬ 
plied  by  the  appropriate  reflection  coefficients.  These 
solutions,  once  developed,  should  prove  useful  for 
providing  building  blocks  to  more  complex  geom¬ 
etries.  For  example.  Figure  25  sketches  some  of  the 
additional  mechanisms  that  might  appear  in  a  ‘fault¬ 
like'  geometry.  The  rays  diffracted  by  the  upper 
corner  appear  after  reflection  from  the  lower  inter¬ 
face  to  be  emanating  from  the  image  shown.  Of 
course,  improved  values  arc  needed  for  the  fields  dif¬ 
fracted  by  the  upper  corner.  Evaluating  the  reflection 
of  the  trapped  wave  also  would  require  'additional 
study.’  These  would  probably  give  fields  of  lower 
value  because  of  loss. 


CONCLUSIONS 

The  scattering  attenuation  function  has  been  de¬ 
fined  for  the  scattering  from  two-dimensional  targets 
using  finite  size  antennas.  Some  results  have  been 
presented  for  the  SAF  of  targets  in  a  lossy  medium. 
Possible  extensions  of  existing  methods  for  this  pur¬ 
pose  have  also  been  discussed. 

The  aulhurs  arc  indchleJ  to  W  D  Burn¬ 
side  for  helpful  suggestions  and  permission  to  use  some  of  hii 
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ment  of  I'.lectrical  Imgineering.  I'lectroScience  Laboratory.  Ohio 
Stale  Unis'ersily.  Columbus.  Ohio,  and  in  part  by  the  Army  Re¬ 
search  Omce.  contract  DAAO29-79.Cm0«2 
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PROPAGATION  OF  SURFACE  WAVES  ON  A  BURIED 


COAXIAL  CABLE  WITH  PERIOOIC  SLOTS* 

J.H.  Richmond 

The  Ohio  State  University  ElectroScience  Laboratory 
Department  of  Electrical  Engineering 
Columbus,  Ohio  43212 

ABSTRACT 

Consider  a  horizontal  coaxial  cable  with  periodic  slots  in  the 
outer  conductor.  This  "leaky  cable"  is  buried  in  the  earth  as  one 
component  of  an  intruder  detection  system.  We  develop  the  theory  for 
surface-wave  propagation  on  the  cable  in  the  presence  of  the  planar- 
air-earth  interface.  Numerical  results  are  included  for  the  phase 
velocity  and  attenuation  constant  as  functions  of  the  various 
parameters.  Data  are  presented  for  the  electric  field  strength  at  the 
air-earth  interface,  and  the  electric  field  distribution  in  the  air 
region  above  the  buried  cable. 
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